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ABSTRACT 

Context. Accretion models predict that fluorescence lines broadened by relativistic effects should arise from reflection of X-ray 
emission onto the inner region of the accretion disc surrounding the central black hole of active galactic nuclei (AGN). The theory 
behind the origin of relativistic lines is well established, and observational evidence from a moderate number of sources seems to 
support the existence of these lines. 

Aims. The aim of this work is to establish the fraction of AGN with relativistic Fe K„ lines, and study possible correlations with 
source physical properties. 

Methods. An XMM-Newton collection of 149 radio-quiet Type 1 AGN has been systematically and uniformly analysed in order to 
search for evidence of a relativistically broadened Fe K„ line. To enable statistical studies, an almost complete, flux-limited subsample 
of 31 sources has been defined by selecting the FERO sources observed by the RXTE all-sky Slew Survey with a count rate in the 3- 
8 keV energy band greater than 1 cts/sec. The 2-10 keV spectra of the FERO sources where compared with a complex model including 
most of the physical components observed in the X-ray spectra of Seyfert galaxies: a power law primary continuum modified by non- 
relativistic Compton reflection and warm absorption, plus a series of narrow Fe line reflection features. 

Results. The observed fraction of sources in the flux-limited sample that show strong evidence of a relativistic Fe K„ line is 36%. 
This number can be interpreted as a lower limit to the fraction of sources that present a relativistic broad Fe K„ line in the wider AGN 
population. The average line equivalent width (EW) is of the order of 100 eV. The outcome of the fit yields an average disc inclination 
angle of 28 ± 5° and an average power-law index of the radial disc emissivity law of 2.4 ± 0.4. The spin value is well constrained only 
in 2 cases (MCG-6-30-15 and MRK 509); in the rest of the cases, whenever a constraint can be placed, it always implies the rejection 
of the static black hole solution. The Fe K^ line EW does not correlate with disc parameters or with system physical properties, such 
as black hole mass, accretion rate, and hard X-ray luminosity. 

Key words, quasars: emission lines - galaxies: nuclei - galaxies: active - X-ray: galaxies - line: profiles 
1. Introduction metal atoms (e.g. George & Fabian 1991*). From a combination 



of fluorescence yield (proportional to the fourth power of the 

The standard scenario for the X-ray emission in active galac- ^^^mic number) and cosmic abundance, the most prominent is 

tic nuclei (AGN) assumes that the observed X-ray power-law ^^^ ^ jj^^ ^^^^^^ ^y ^^^j^^j ^^^ ^^ ^^ j^^y. 
continuum originates in the inner regions of the AGN closest 

to the central super-massive back hole, via inverse Compton When the Compton reflection originates in distant mate- 
scattering of soft-energy photons in a corona of relativistic rial Uke the mol ecular torus envisaged in unification models 
electrons located somew here above the accretion disc (e.g. jAntonucci 19931 ). the profile of the Fe K„ emission fine is nar- 
IHaardt & Maraschi 1993l l. The X-ray illumination of optically row and unresolved by present X-ray detectors. In contrast, the 
thick cold matter, such as the molecular torus and/or the accre- reflection component originating in the inner accretion disc is 
tion disc, gives rise to a Compton reflection spectral component affected by the black hole's strong gravitational field, which 
containing a series of fluorescent lines from K^ transitions in modifies the line profile. The numerous studies of the effect 

of gravity on the narrow emission line all agree by describing 

Send offprint requests to: I. de la Calle Perez the resulting relativistic profile as skewed and asymmetric be- 
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cause of a combination of kinematic and relativistic effects (e.g. 
IFabian et al. 2000l l. with a red wing extending towards low X- 
ray energies due to gravitational redshift and transverse Doppler 
redshift (e.g. lFabianet al. 1989l l. 

Since relativistic lines originate within a few gravitational 
radii from the central object, the study of their shape and inten- 
sity may represent a potential probe of the physical processes 
taking place in the innermost regions of the AGN. The line pro- 
file is in fact very sensitive to the accretion disc properties, such 
as the radial extension and dependence of the line emissivity, the 
ionisation state of the material, the observer's inclination angle 
relative to the disc axis, and the spin of the black hole (for a re- 
view see Reynolds & Nowak 2003 IFabian & Miniutti 2005 1 and 
IMatt 2006l and references therein). 

The first unambiguous observational evidence of relativis- 
tic Fe Kq, lines in the X-ray spectra of AGN was found in 
the spectrum of MCG- 6-30-15 obtained by the ASCA satel- 
lite dTan aka et al. 1995") and lately, in several other Seyfert 
Galaxies (Nandraetal. 1997) to the point where the broad 
line was considered a common feature of AGN. Nevertheless, 
the presence of relativistic broad Fe Ko, lines is nowadays 
more controversial than ever since only a handful of sources 
seem to possess a truly relativistically broadened Fe K„ 
line, e.g. MCG-6- 30-15 (Fabian et al. 2002 ) or NGC 3516 
jMarkowitz et al. 20 08). When studies based on sizable 
samples of sources are considered, the average fraction of 
sources with relativistic broad Fe Kq, lines is never higher 
than 40% (Porquet et al. 2004[ IJimenez-Bailon et al. 20051 



IGuainazzi et al. 2006, Nandra et al. 20071 1. 

In the past decade the advent of Chandra, XMM-Newton 
and Suzaku has undoubtedly provided deeper knowledge and 
more thorough understanding of AGN X-ray spectral properties 
by means of superior spectral resolution and higher through- 
put. For instance, it has been confirmed that about half of 
quasars and Seyfert galaxies show the so-called warm ab- 
sorbers, i.e., ionised material outflowing along the line of sight, 
which is revealed as a series of absorption features imprinted 
mainly in the soft X-ray spectral band ( |Piconcelli et al. (2005)[ 
IBlustin et al. 20051 and references therein). Emission lines from 
highly ionised atoms, including iron, can also be observed (see 
IBianchi et al. 2005 ). The presence of this absorbing/emitting gas 
may introduce further complexities into the Fe K band: it has 
been shown that warm absorbers with sufficiently high col- 
umn density and ionisation state may distort the continuum 
underlying the iron line, mimicking a relativistic broad red 
wing (e.g. NGC 3783, Reeves et al. 200 4j and M CG-6-30-15, 
IMilleretal. 20081 see also ,Turner & Miller 20091 for a review). 
Also, a blend of emission lines from helium and hydrogen-like 
iron at 6.7-6.97 keV can be mistaken for a broad disk line ob- 



served at high inclination (NGC 7314, Yaqoob et al. 2003 and 
MRK 590, [Longinotti et al. 2007[ ). While the latter scenario is 
fully taken into account in our baseline model, only a simple pa- 
rameterization of the warm absorber was adopted, because an 
unknown number of a priori absorbing systems do not allow a 
systematic analysis to be performed (see Section [33] for details). 
The FERO (finding extreme relativistic objects) project is 
part of a wider investigation on AGN carried out on archival 
XMM-Newton data (see Section |2|. It was designed to address 
two fundamental questions on the relativistic broad Fe K^ fine: 

i) how common is relativistic broadening in AGN, and 
ii) does the presence of a broad line depend on other source's 
physical properties ? 



To try to answer these questions, all the sources in our sam- 
ple were fitted with one and the same model, which includes 
all absorption and emission components known to be potentially 
present in AGN (see Section [33] ) and able to affect the emission 
in the iron K band. 

This paper reports the results of the spectral analysis on 
the individual sources of the sample. A companion paper 
(Longinotti et al. in preparation) is devoted to the analysis 
of stacked spectra. The structure of this paper is as follows. 
Section |2] describes the selection of the sample. Section [3] de- 
scribes the analysis procedure, including the spectral analysis 
and the spectral model used to describe our data. Section|4]sum- 
marizes our results to be discussed in Section|5] Section|6]closes 
the paper with the conclusions drawn from the study presented 
here. 

2. The Sample 

The FERO AGN collection proceeds from the CAIXA catalogue 
of AGN recently published by Bianchi et al. (2009a,b). CAIXA 
is the largest catalogue of high signal-to-noise X-ray spectra of 
AGN which consists of all the radio-quiet X-ray unobscured 
(A^H < 2 X lO^^cm"^) AGN observed by XMM-Newton in tar- 
geted observations, whose data are public as of March 2007. The 
sample, through a complete and homogeneous spectral analy- 
sis, is characterized in terms of the parameters adopted by the 
best-fit models. CAIXA includes a total of 77 quasars and 79 
Seyfert galaxies. The redshift distribution spans from z=0.002 
to z=4.520 (almost 90% within z<l) and the distribution of the 
hard X-ray luminosities covers a range between L2-io/t«v = 2.0 
10**' - 3.9 10'*'' erg s"'. Since the source selection criteria used 
in this work are by large the same described there, the reader 
is deferred to these works for more details. Here, those general 
aspects that are more relevant to the relativistic Une analysis car- 
ried out in this work are recalled. 



2.1. The FERO AGN Collection 

The starting sample consisted of 161 radio-quiet X-ray unob- 
scured sources targeted by XMM-Newton with public data up 
to April 2008. Only sources with local column density from 



cold gas lower than Nh < 2 x 10 cm are included in the 
sample to avoid heavily absorbed spectra in the 2 to 10 keV 
spectral region. With respect to CAIXA (which sums up 156 
sources), the following 5 AGNs were added: ESO 511-G030, 
IRAS 05078-H1626, MRK 704, NGC 3227, NGC 526A. The 
first three sources were specifically proposed and granted in the 
XMM-Newton A06 for the purpose of the FERO project; the 
other two objects became publicly available in the data archive 
during the development of the FERO analysis. NGC 4151 full- 
fills the required conditions to be in our sample, but it was ex- 
cluded (both from FERO and CAIXA) based on its complex 
X-ray spectra and extreme spectral variability attributed to sev- 
eral absorbing systems dPuccetti et al. 2007 ). Based on a signal- 
to-noise criterion, as reasoned in Section 13.21 12 sources are 
dropped from the starting sample, leaving what constitutes the 
FERO sample with a total of 149 sources. 

Out of the 149 sources, 67 are classified as quasars (RQQs) 
and 82 classified as Seyfert Is (Sy) according to the value 
of the absolute optical magnitude M^ (QSO Mb < -23; Sy 
Mb > -23) as defined in the Veron-Cetty & Veron catalogue 
2006 ( Veron-Cetty & Veron 2006) . Radio Loud objects are ex- 
cluded from the sample according to the value of the radio- 
loudness parameter (J?; IStocke et al. 19921 ). For QSO the condi- 
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tion log(R) > 1 is applied to reject a source, while in the case of 
Seyferts, in addition to log(R) > 2.4, sources are also excluded 
if log(Rx) > -2.755 (Panessaet al. 2007), where R x is th e X- 
ray radio-loudness parameter (JTerashima & Wilso n 20031) (see 
IBianchi et al. 2009a| for details). No redshift restriction has been 
imposed in our sample, which contains mostly local sources with 
90% of the sources at a redshift < 0.5 and 60% at a redshift 
< 0.1. A distinction has been made between Broad Line (BL; 
40% of full sample) and Narrow Line (NL; 20% of full sam- 
ple) sources, using as threshold the value of the Full Width Half 
Maximum (FWHM) of the H^g line whenever available (65% of 
the sources in our sam ple): H^g > 2000 km/sec for BL and H^ < 
2000 km/sec for NL (Osterbrock & Pogge 198 7). This distinc- 
tion between NL and BL refers only to this standard limit on 
H^, while no optical Type 2 objects are present in the sample. 
When needed, and not derived in this work, source properties 
are extracted from the CAIXA catalogue. 



2.2. The Flux-limited Sample 

The selection criteria with which CAIXA and FERO were as- 
sembled do not provide a complete sample in the sense that the 
sources are not selected according to a physical property: the 
main point of the FERO project is to include as many objects 
as possible and this of course depends predominantly on the 
availability of public observations. Because of the nature of the 
FERO project, one of the fundamental requirements was to iden- 
tify a complete, unbiased subsample of sources with high signal 
to noise in order to derive meaningful constraints on the prop- 
erties of the (unknown) parent population of local radio-quiet 
AGN. Hence, sources from the RXTE all-sky Slew Survey (XSS, 
[Revnivtsev et al. 2004 ) having a count rate in the 3-8 keV energy 
band greater than 1 cts/sec and fulfilling the FERO source selec- 
tion criteria were identified. This defines a flux-limited sample of 
33 sources. The XSS is nearly 80% complete at the selected flux 
level for sources with Galactic latitude greater than 10°. For two 
of them, UGC 10683 and ESO 0141-G055, no XMM-Newton 
data were available as of April 2008. This leaves the number of 
XSS-selected bright sources in the FERO sample to 31 (listed 
in Table IaTI in Appendix lAli. Throughout the paper, these 31 
sources will be referred to as the flux-limited sample. 



3. Analysis 

The data corresponding to the 149 observations have been 
uniformly analysed using SASv6.5 JGabriel et al. 2004) and 
the latest calibration files available. Event lists were obtained 
for the EPIC-pn camera following standard SAS data reduc- 
tion procedures. Due to its larger effective area, only EPIC-pn 
dStriider et al. 20011 1 data has been considered in this work. The 
filtering of event lists for periods of high background activity 
was performed by maximizing the source signal to noise as in 
|Piconcelli et al. (2004)| The typical source extraction regions are 
circular and of the order of 50"in radii for small window mode, 
while for full frame and large window modes, the source extrac- 
tion regions range from 20 to 40"in radii. Background spectra 
were extracted from circular source-free regions, except for ob- 
servations taken in small window mode where blank-field event 
lists were used as described in [Read & P onman ( 2003)| of a 
50"in radii and close to the target. Spectra affected by a pileup 
larger than 1 % were rejected. Full details of the data reduction 
and spectra accumulation are reported in Bianchi et al. (2009a). 



3.1. XMM-Newton observations 

The data of the observations presented here were public as of 
April 2008. All the observations are target observations with ex- 
posure times ranging between 1 ksec and 400 ksec, with 90% of 
the observations below 100 ksec. 



3.1.1. Multiple observations 

When multiple observations of the same source were available, 
the individual spectra were combined and treated, for all pur- 
poses, as a single observation. The FTOOLS tasks mathpha, 
adnnf, and addarfwere used correspondingly to sum the source 
and background counts spectra and the response matrix and an- 
cillary files, using the exposure times of the individual observa- 
tions to be summed as weights in these last two cases. However, 
spectra were only combined if the observations were taken with 
the same observing mode and if the source was in similar flux 
state, i.e., if the total flux and power-law spectral index in the 
2-10 keV energy band were consistent within the statistical er- 
rors. If these criteria were not met, then the observation with the 
longest exposure time was chosen to prevent introducing source- 
related bias in the selection. This procedure ensures that mixing 
different spectral states from a given source does not affect the 
results derived from our data set. 

A total of 22 sources in our sample have multiple observa- 
tions. Table [1] lists on a source-by-source basis XMM-Newton 
observation IDs that have been combined. For completeness. 
Table IbTI in Appendix IBJ lists those sources for which multiple 
observations are available, but only one was considered; so no 
spectra were combined. This approach is different from the one 
taken in CAIXA, where for sources with multiple observations, 
the one with the longest exposure was selected systematically as 
independent of the flux state. The importance of achieving high 
signal-to-noise spectra for the FERO project, meant that it was 
preferred to co-add multiple observations of the same source, 
when available. 



3.2. Spectral analysis 

The time-averaged spectra were re-binned in order not to over- 
sample the intrinsic energy resolution of the EPIC-pn camera 
(^150 eV at 6 keV) by a factor larger than 3, while making 
sure that each spectral channel contains at least 25 background- 
subtracted counts to ensure the applicability of the;i'^ goodness- 
of-fit test. Fits were performed in the 2-10 keV energy range. 
To ensure good quality spectral fits, only those source spectra 
with at least 17 d.o.f. in this energy range were considered in the 
analysis. (12 out of the 161 sources in the initial sample did not 
meet this criterion.) For the spectral analysis and fitting, XSPEC 
vl2.3.0 (Arnaud 1996) was used. Throughout the analysis, solar 
abundances were assumed after Anders & Ebihara (1982). 



3.3. Spectral model 

All spectra were fitted with the following baseline model: 



e-A'n . w(r, Ni^, ^) ■ A[£-r + c(r, R) + Y^ Gi + Ky(9,f5, a)] (1) 

(=1 

where the diff'erent model components are 
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Table 1. List of sources within the FERO sample where multiple observations are available. 



Source 


List of 


Observation 


Observation 


Source 


List of 


Observation 


Observation 


Name 


XMM-Newton 


Date 


Exposure 


Name 


XMM-Newton 


Date 


Exposure 




Observation 




Time 




Observation 




Time 




IDs combined 




(ksec) 




IDs combined 




(ksec) 


PG1440-H356 


0005010101 


2003-01-01 


17.2 


1H04 19-577 


0148000401 


2003-03-30 


11.1 




0005010201 


2003-01-04 


10.6 




0148000501 


2003-06-25 


10.7 




0005010301 


2003-01-07 


18.1 




0148000601 


2003-09-16 


11.3 




0107660201 


2001-12-23 


24.2 




0148000201 


2002-09-25 


11.8 


AKN564 


0006810101 


2000-06-17 


7.4 




0148000301 


2002-12-27 


7.6 




0006810301 


2001-06-09 


7.5 




0148000701 


2003-11-15 


1.1 




0206400101 


2005-01-05 


69.2 




0112600401 


2000-12-04 


5.7 


MCG-6-30-15 


0029740101 


2001-07-31 


55.8 


PGl 115-1-080 


0203560201 


2004-06-10 


65.5 




0029740701 


2001-08-02 


85.7 




0203560401 


2004-06-26 


71.2 




0029740801 


2001-08-04 


86.8 




0082340101 


2001-11-25 


53.8 




0111570101 


2000-07-11 


28.8 


IRAS 17020+4544 


0206860101 


2004-08-30 


13.4 




0111570201 


2000-07-11 


37.9 




0206860201 


2004-09-05 


12.3 


NGC5548 


0089960301 


2001-07-09 


58.7 


MRK507 


0300910401 


2005-06-17 


15.8 




0109960101 


2000-12-24 


16.0 




0300910701 


2005-09-01 


11.4 




0089960401 


2001-07-12 


19.8 


MRK279 


0302480401 


2005-11-15 


41.4 


MRK876 


0102040601 


2001-04-13 


2.6 




0302480501 


2005-11-17 


40.7 




0102041301 


2001-08-29 


2.4 




0302480601 


2005-11-19 


22.8 


ESO15-1G011 


0103861701 


2000-09-29 


5.0 


MRK766 


0109141301 


2001-05-20 


89.6 




0103862001 


2001-10-31 


4.5 




0304030301 


2005-05-25 


68.9 


MCG-01-13-025 


0103863001 


2002-08-28 


4.3 




0304030401 


2005-05-27 


65.8 




0103861401 


2000-08-30 


1.5 




0304030501 


2005-05-29 


64.9 


NGC3516 


0107460601 


2001-04-10 


64.3 




0304030601 


2005-05-31 


63.1 




0107460701 


2001-11-09 


88.7 




0304030701 


2005-06-03 


20.4 


REJ2248-511 


0109070401 


2000-10-26 


10.1 




0096020101 


2000-05-20 


25.7 




0109070601 


2001-10-31 


9.8 




0304030101 


2005-05-23 


66.2 


H0557-385 


0109130501 


2002-04-03 


4.0 


SDSSJ135724.51+652505.9 


0305920301 


2005-04-04 


21.1 




0109131001 


2002-09-17 


6.5 




0305920601 


2005-06-23 


12.0 


TONS 180 


0110890401 


2000-12-14 


20.3 


LBQS1228+1116 


0306630101 


2005-12-13 


60.4 




0110890701 


2002-06-30 


12.6 




0306630201 


2005-12-17 


83.2 


NGC7469 


0112170101 


2000-12-26 


12.3 


MRK509 


0306090201 


2005-10-18 


59.8 




0112170301 


2000-12-26 


16.1 




0306090301 


2005-10-20 


32.4 




0207090101 


2004-11-30 


59.2 




0306090401 


2006-04-25 


48.6 




0207090201 


2004-12-03 


55.0 




0130720201 


2001-04-20 


27.8 


NGC3783 


0112210101 
0112210201 
0112210501 


2000-12-28 
2001-12-17 
2001-12-19 


26.1 
51.2 
93.2 




0130720101 


2000-10-25 


20.7 



Notes. Multiple observations were combined where available and treated, for all purposes, as a single observation. Observation IDs in bold 
correspond to those observations that were combined, while observation IDs not marked in bold were discarded. 



Galactic absorption (e"°"'^"), where A^h was fixed 
to the galactic column density and cr is the pho- 
toelectric cross section of the process according to 
[Morrison & McCammon (1983)t 

Power law (N(E) oc E~' ), where F is the photon index of the 
primary power-law X-ray source spectrum; 
Compton reflection from neutral material (C(F, R)), where R 
is the reflection factor as compared to material emitting over 
In in solid angle (see however Murphy & Yaqoob 2009 for 
some caveats about this interpretation). The high-energy cut- 
off of the primary spectrum was fixed to 100 keV. The view- 
ing angle of the reflecting material was fixed to 18°, while 
the reflection normalization (R) was left free during the fit- 
ting procedure. The model used was the XSPEC PEXRAV 
component (Magdziarz & Zdziarski 1995|l; 
Ionised absorption {W{r,N'^^,^)) from warm gas in the lo- 
cal AGN environment was introduced through the XSPEC 
model ABSORl. Although more accurate and complex mod- 
els exist for warm absorption, they would overfit the lower- 
quality data in our sample when uniformly applied to the 



whole sample. Both the intrinsic column density (A^^) and 
ionisation parameter (^) were left free during the fitting pro- 
cedure, while the continuum slope (F) was tied to the photon 
index of the primary power-law continuum. The temperature 
of the absorber was fixed to 2x10^ K (~17 eV), represen- 
tative of the temperature found in warm absorbers (see e.g. 



Krongold et al. 2007) ; 



Narrow lines (Gi): 4 zero-width Gaussian fines were in- 
cluded with centroid energies fixed at 6.4 keV (neutral 
Fe I K„), 6.7 keV (ionised Fe XXV), 6.96 keV (ionised 
Fe XXVI), and 7.06 keV (neutral Fe I K^). The normal- 
ization of the Fe I K^ line was hmited to less than 16% 
that of flie Fe I K^ line flux jMolendi et al. 20031 1. The 
two ionised lines were included in the baseline model 
as several studies show that they represent an important 
spectral component in AGN. Among the possible phys- 
ical interpretations of these features are an origin in the 
accretion disc (e.g. IReeves et al. 20011 [Pounds et al. 20011 
and IMiUer et al. 20071) or produced in photoionised 
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circumnuclear matter (e.g. IBianchi & Matt 20021 and 
IBianchi et al. 2005l l: 

• Fe I 6.4 keV Compton shoulder (Gi). The Fe I K(j, Compton 
shoulder of the narrow component of the Fe I fluorescent 
line was considered by using a Gaussian line at 6.3 keV, with 
0.05 keV as the fixed width and free normalization allowed 
to be no more than ~20% that of the 6.4 keV Fe I Kq, line, 
following Matt (2002); 

• Relativistic line {Ky(6,/3,a), Kyrline model). To model the 
effect of the strong gravitational field on the emission re- 
gions of the disc closest to the central black hole, the Ky 
set of models within XSPEC was used (Dovcia k et al. 2004b . 
The model considers accretion disc emission around a rotat- 
ing black hole with a given spin, a. In contrast to models 
such as diskline or laor, Kyrline allows the spin to be fit- 
ted as a free parameter For the line emission, a fixed cen- 
troid energy of 6.4 keV (in the rest frame of the source), is 
assumed, while the emission between the innermost stable 
orbit and 400 gravitational radii is integrated, keeping these 
limits fixed during the fitting procedure. The radial depen- 
dence of the disc emission is modelled with a power law 
of index fi (oc r"^), while Laor's limb darkening law was 
adopted to characterize the angular dependence. Here, the 
radial disc emissivity is referred to as /3. The disc inclina- 
tion angle (0), yS, and a are free parameters during the fitting 
procedure. 

In summary, the model that was uniformly fitted to our data 
consists of 9 components with a total of 14 free independently- 
fitted parameters. For those cases in which the value of the 
FWHM of the H^ line is available, the model has 15 free param- 
eters (see Section [3.3. II ). The fitting procedure was carried out 
over different steps. First, the power law with neutral reflection 
and absorption were fitted to the data, followed by the addition 
of the narrow line components one by one (in the order listed 
above). The last model component added to the fit was the rel- 
ativistic line. After each addition of a new component, a fit was 
performed, and no model component discarded regardless of the 
goodness-of-fit. This implies that the x^ values obtained from 
the fits correspond to the best-fit given by this model. Its impor- 
tant to stress that no attempt has been made as a whole or on 
a source-by-source basis to find the best-fit model yielding the 
lowest possible ;t'^ value. For reference, table lC.ll gives. for those 
sources belonging to the flux-limited sample, the most relevant 
parameters of the best fit for the model considered. 

3.3.1. Limitations to model parameters 

During the fitting procedure, some limitations are introduced for 
the following parameters of the baseline model: 

o The disc inclination angle in the Kyrline model was limited 
to a maximum value of 60°. In the Seyfert unification sce- 
nario dAntonucci 19931 1. Type 1 objects are not expected to 
have large viewing angles. It has been found that relaxing 
the restriction on the inclination angle yields a number of 
detections of broad and highly skewed lines with equivalent 
widths of the order of 1 keV. This effect is simply because 
in sources with poor statistics above 7 keV, the broad line 
model fits the continuum, also reflected by required large 
disc-inclination angles (>80°). 

o The disc emissivity in the Kyrline model was limited to a 
maximum value of 6 given that steeper profiles are not ex- 
pected, even taking strong general relativistic effects into ac- 



count on the primary emission such as light bending (see 
IMiniutti & Fabian 2004). 
o To possibly constrain the parameter space, limitations to the 
spin parameter in the Kyrline model were tested. Three fit- 
ting runs were produced with black hole spin fixed to 0.998 
{Kerr), fixed to {Schwarzschild), and free to vary between 

and 0.998. In the first two cases, no systematic difference 
was found when comparing the results in terms of x' im- 
provement, leading us to conclude that none of the two ex- 
treme spin values can be assumed a priori. Therefore, the 
results presented are extracted from the run where the black 
hole spin is left free to vary during the fitting procedure. 

o The power-law index (F) and the reflection fraction (R) were 
limited following the work by |Dadina (2008)| on a BeppoSAX 
sample of Type 1 Seyfert galaxies. Given the distributions of 
r and R derived by Dadina (2008) a limit of three times the 
Icr standard deviation of the average values is chosen for F 
and R (F e [0.45:3.33] ; R e [0.:3.33]). 

o When available in the literature, the value of the FWHM 
of the optical Hp line was used to set an upper limit to the 
width of the narrow component of the 6.4 keV neutral Fe 

1 K„ line, leaving the width of the line free. This choice is 
consistent with the possibility that for some sources, the Fe 
I K„ line may be produced in the optical broad line region 
(see e.g. Bianc hi et al. 20081 and references therein). When 
the H/j FWHM was not available, the Fe I K„ line width was 
assumed to be unresolved (cr = 0). 

3.3.2. Test of spectral model components 

During the process of selecting the baseline model, several tests 
were performed with different model components. The most im- 
portant ones are summarised below. 

o To contemplate the hypothesis that the reprocessed AGN 
continuum in Type 1 objects is made by contribution of both 
torus and disc reflection, the presence of two reflection com- 
ponents at the same time was tested. The second PEXRAV 
disc component was relativistically blurred using the convo- 
lution model kdblur, to take possible relativistic effects into 
account. After close inspection of the contribution from each 
physical system, it was concluded that the sensitivity band 
of XMM-Newton does not allow disentangling the reflection 
component from the torus from that of the disc. Since the 
presence of torus reflection is very well supported obser- 
vationally by the almost ubiquitous narrow Fe I K„ com- 
ponent (IBianchi et al. 2004l[Yaqoob & Padmanabhan 2004| l, 
only this PEXRAV component has been included in our base- 
line model. 

o The effect of reflection from an ionised accretion disc was 
tested by replacing the torus reflection with the model com- 
ponent PEXRIV ([Magdziarz & Zdziarski 1995[): Eq.[T]is thus 



modified to include the disc ionisation parameter ^ in the 
function C(r, R), which is left free during the fitting proce- 
dure. The centroid energy in the Kyrline model was consis- 
tently set to 6.7 keV. The results from this test are discussed 
in Sects. |4]and|5] 

4. Results 

The EW of the Kyrline component was used as a proxy for estab- 
lishing significant detections of relativistic Fe I Kq, lines in the 
sample. The line EW is defined in this work as the ratio of the 
flux of the Kyrline component to the flux of the continuum, as 
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defined in Section [331 integrated over the energy range between 
±5% around the fine rest-frame energy (6.4 keV). The error in 
the normalization of the Kyrline component is used to estimate 
the uncertainty in the corresponding line EW. The threshold for 
significant detections is set to 5cr confidence level (c.l.). This 
is a somewhat more conservative approach than normally taken 
when reporting relativistic line detections in individual objects, 
but it has to be kept in mind that the FERO project is a statisti- 
cal study of a large sample that includes spectra of very different 
levels of statistical quality. For those sources with no significant 
line detection, the upper limits to the line EW at the 90% c.l. are 
provided. In the following, errors, upper and lower limits, on the 
relevant parameters are given at the 90% c.l. for one interesting 
parameter (Ax^ = 2.71). 

4. 1 . Sources with detected broad Fe lines. 

Figure [1] shows the EW of the broad Fe Kq. line at 6.4 keV as 
a function of source counts in the 2-10 keV band. This plot 
shows two things. First, significant broad Fe K„ line detections 
are concentrated in the portion of the plot dominated by spectra 
with good statistical quality (roughly with > 1.5x10^ hard X-ray 
counts). Second, overall, the distribution of line EW upper limits 
versus hard X-ray counts follows the expected trend; i.e., in the 
absence of a significant broad line, the sensitivity (the ability to 
detect a broad line) goes as (l/V(cf52-i{)kev)- The dispersion of 
EW upper limits around this direction, for a given range of hard 
X-ray counts, could be explained in terms of scattering of the 
model parameters. 

Table |2] provides the detailed list of the significant relativis- 
tic broad Fe K„ line detections with the best-fit parameters of 
the Kyrline model. From now on, this table is referred to as the 
list of detected relativistic Fe K„ lines in the FERO sample. In 
the fitting procedure, relativistic effects from neutral and ionised 
reflection were tested as explained in Section [3.3.2l The test for 
a neutral Fe K^, broad line yields in total 1 1 significant detec- 
tions at >5cr C.L. The fit with the model with an ionised Fe Kq, 
line yields a detection in all the 1 1 sources, except IC4329A, 
1H0707-495, and MCG-5-23-16, which only provide detection 
in the neutral case. In the 8 sources where a relativistic broad 
Fe Kq, line is equally detected by means of neutral and ionised 
reflection models, the best-fit models were carefully examined 
and the model with the lowest x^ preferred. MRK 509 is the 
only source for which the ionised reflection provides a signifi- 
cant better-fit statistic (with the probability of the improvement 
being by chance of 2.2x10"^). MRK 766 and ARK 120 on the 
contrary, show significant broad Fe lines only in the ionised case. 
At the end, neutral and ionised runs were merged, which yields 
13 sources with strong evidence of a broad Fe line. Last, it is 
worth pointing out that the only source where a relativistic broad 
Fe K„ line is claimed for the first time in this work is ES051 1- 
G030, with an EW of 57+f/ eV. 

Two sources, PG121 1 + 143 and 1H0707-495, are charac- 
terised by extremely high EW of the Fe Kq line (Figure [1] see 
also Table |2]l. The data included in our sample for both sources 
have similar spectral shapes, showing deep spectral features at 
~7 keV and spectral curvature in the Fe energy band. Pounds 
et al. (2003a) argue that PG121 1 + 143 presents a high column 
density (5 x 10^^ cm"^) of highly ionised matter along the line 
of sight partially covering the central hard X-ray source and 
that the apparent presence of a relativistic Fe Kq emission line 
could be an artifact of absorption. However, the relativistic broad 
line interpretation is not discarded, and the same authors find 
that a broad line with a large EW, of the order of 600-900 eV, 



and a combination of absorption features, provide a good in- 
terpretation of the data. This value of the Fe Kq line EW is in 
good agreement with the value found in this work (635^24? ^^)- 
Nevertheless, whether there is a relativistically broadened Fe Kq 
line in this source is debatable, as new results from a more re- 
cent analysis on data from 2007 indicate (IPounds et al. 2009] l. 
1H0707-495 also shows a prominent flux drop around 7 keV, 
which can be modelled with an absorption edge assuming large 
Fe overabundance (Gal lo et al. 20041 1. Fabian et al. (2004) inter- 
pret the drop around 7 keV in terms of relativistically blurred 
ionised reflection from the accretion disc, and derive an EW for 
the broad Fe Kq line of 1800 eV, in agreement with the value 
found in the FERO analysis (1775!^'^ eV). Fabian et al. (2004) 
proposed that the X-ray spectrum of this source is disc reflection 
dominated, and as the high value of the line EW implies, high el- 
ement abundances were required, in a similar way to absorption- 
dominated models (IGallo et al. 2 004). Recently, Fabian et al. 
(2009) have presented new timing-based arguments in favour of 
the relativistic nature of the broad features in this source with the 
simultaneous detection of broad Fe L and K lineo However, 
Miller et al. (2010) has recently claimed that the Fabian et al. 
(2009) interpretation of the timing properties of the source is not 
unique and that solutions that do not invoke relativistic repro- 
cessing might still be viable. 

Figure |2] shows the ratios of the data to the continuum and 
data to the best-fit model in the 2-10 keV energy range for 
the sources with a relativistic 6.4 keV Fe Kq broad line detec- 
tion. For MRK509, MRK766 and ARK120, the figures were 
produced with the model including the ionised continuum and 
6.7 keV Fe Kq line. The latter set of plots indicates how well the 
baseline model used in this work describes the data. The devi- 
ation of the best-fit model and the data is limited below ~10% 
in the whole 2-10 keV energy range used for fitting, and below 
~5% if one excludes the last energy bins (>8 keV) with poorer 
statistics. Despite the fact that by selection the sources in FERO 
are unobscured (A^h < 2 x lO^^cm"^) for two sources, MCG-5- 
23-16 and IC4329A, a downturn in the ratio plot below 3 keV 
is noticeable and could indicate extra cold absorption, which in 
principle would not have been taken into account by our baseline 
model. However, for the particular case of these two sources, it is 
believed that this downturn cannot be attributed to extra cold ab- 
sorption since the parameters derived from the warm absorption 
model (ABSORI) used yield values compatible with those ex- 
tracted from CAIXA (as seen in TableEU, A^h (z=0) (1.20!°-^^) 
X 10^2 cm-2, and (0.33!°°^) ^ 10^2 cm^^ for MCG-5-23-16 and 
IC4329A, respectively, and in both cases, compatible with neu- 
tral absorption (^ < 0.02 erg cm"' s"' for MCG-5-23-16 and 
^ < 0.47 erg cm"' s"' for IC4329A). For these two sources, any 
extra neutral absorption that is present should have been taken 
into account by the model component ABSORI. Table |2] lists in 
the last column the x^ and corresponding number of degrees of 
freedom for those sources with evidence of a relativistic Fe Kq 
line. 



4.2. Fraction of relativistic Fe lines 

4.2.1. Detection fraction in FERO 

The fraction of relativistic Fe lines detected in the FERO sam- 
ple is 9% (13/149). Considering only the sources in the flux- 



' We point out that the XMM -Newton spectra of 1H0707-495 in- 
cluded in the FERO sample correspond to the ones used by Fabian et al. 
(2004) and that Fabian et al. (2009) used a different dataset. 
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Fig. 1. Equivalent width of the relativistic broad 6.4 keV Fe Kq, line vs. hard X-ray counts (2-10 keV) for the FERO sources. Filled 
circles indicate line detections at >5cr confidence level (where error bars indicate the 90% confidence level intervals), while filled 
squares indicate line upper limits at the 90% confidence level. White stars indicate sources belonging to the flux-limited sample (see 
Section l2?2l for details). 



limited sample, the detection fraction rises to 36% (11/31). 
These two numbers do not provide specific information on the 
intrinsic fraction of AGN that have broad Fe lines. They are lim- 
ited by FERO being made of spectra of disparate quality and 
by the unavailability of a well-defined complete AGN sample. 
Nevertheless, the observed detection fraction can be considered 
as a lower limit for the intrinsic number of AGN that would show 
a broad Fe line if, for example, all sources were observed with 
the same signal-to-noise. 

It is possible to estimate statistical errors on the detection 
fractions to reflect the uncertainty introduced by considering a 
limited number of sources. Detection fractions are calculated as 
vali/vah, where val\ and vah indicate a given number of sources 
used to work out a given detection fraction. For the errors of 
vail 2, Poisson statistics are assumed, where the errors are calcu- 
lated as the average of the upper (Aul) and lower limit (All) of 
the 84% c.l. (which correspond to the commonly used Aul = 1 + 
'\Jvalj2 + 0.75 and All = -^Jvaljj - 0.25, according to Gehrels 
(1986)), both of which are good approximations for this moder- 



ate confidence level and are accurate to better than 5% even for 
low values of vali^ (see Gehrels 1986 for more information). 
The error of val\/val2 is propagated quadratically. In this way, 
the detections fractions and their associated errors are 9±3% and 
36±14% for the FERO sample and for the flux-limited sample, 
respectively. 

A broader understanding of the incidence of the relativis- 
tic line in our work can be achieved by considering the upper 
limits on the line intensity that characterize the majority of the 
sources. Within the FERO sample, ~16±4% of broad line EW 
upper limits are below 100 eV. This value corresponds to the 
EW of a neutral disc with an inclination of 60°, for an illuminat- 
ing power law of r=2 and solar abundances (IMatt et al. I992l l. 
If the analysis is restricted to the 20 sources with no significant 
broad line within the flux-limited sample, the fraction of sources 
in which the fine EW is limited below 100 eV is 50% (10/20) 
(for reference, for 75, 50 and 25 eV, the fractions are respec- 
tively 25%, 20%, and <13%). However, 5 sources within the 
flux-limited sample have measured EWs lower than and incon- 
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Fig. 2. Set of two figures for those sources with a relativistic 6.4 ke V (6.7 ke V) Fe Ka hne EW detection >5cr confidence level: data 
to the best-fit continuum ratio (left set column), and data to the best-fit model ratio (right set column). The vertical dashed line is 
placed for reference at 6.4 keV (6.7 keV) in the rest frame of the source. On top of each figure the model used for the ratio, the 
corresponding source name, and the observation ID are given. For those sources with summed observations a dummy observation 
ID is given. 
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Table 2. Sources within the FERO sample where a relativistic Fe Kq, line is detected with a significance >5cr. 



Source 


Type*'* 


L^-IO kcV 

(10«ergs-i) 


Qj,2-IOkeV 

(10^ cts) 


EW 

(eV) 


d 

n 


a 


P 


;r7dof 


References*^* 


IC4329A 


BLSY 


56.4 


9.123±0.009 


39+'* 

36+" 

57+^" 
J/ [J 


28+°, 


>0.0 


<1.3 


238.1/161 


(1), (2), (3), (4) 


MCG-5-23-16 


NCSY 


14.2 


8.853±0.009 


Zl 3 


>0.0 


<1.6 


271.4/162 


(5), (6) 


ESO511-G030 


NCSY 


22.4 


1.740±0.004 


18+:^ 


>0.0 


<1.1 


198.2/161 




MCG-6-30-15 


BLSY 


5.4 


13.927±0.012 


260+1' 
260+=° 

997-f4l 
^"^'-78 


40!^ 


oMroi 


4 j+0.2 

2 9+''-^ 

-M 

2 8+0-2 

Z.0_Q3 


227.6/161 


(7), (8), (9), (10) 


NGC4051 


NLSY 


0.3 


2.040±0.004 


22+'' 


>0.46 


198.8/162 


(11), (12), (13) 


NGC3516 


BLSY 


3.6 


2.810±0.005 


27!^ 


>0.48 


223.1/161 


(14), (15), (16) 


NGC3783 


BLSY 


11.4 


9.440±0.010 


88!1] 


<8 


>0.16 


9 7+0.1 
A.I ^2 


304.5/161 


(17), (18) 


NGC3227 


BLSY 


1.2 


3.541±0.006 


66!l' 


23+^ 


>0.0 


1 Q+0-6 
05 


230.9/162 


(19), (20) 


MRK509 (*) 


BLQ 


103.5 


7.471±0.009 


170!?^ 


53!| 


0.78t°l 


>3.8 


277.3/158 


(21), (22) 


MRK766 (*) 


NLSY 


6.1 


7.816±0.008 


234^^3 


20+3 


>0.47 


2 7+0-2 

^■'-0.1 


300.3/160 


(23), (24) 


ARK120 (*) 


BLSY 


90.5 


3.473±0.006 


140!3[ 


>59 


>0.0 


9 9+0.6 


212.3/158 


(25) 


1H0707-495 


NLSY 


46.8 


0.110±0.001 


1775!^^i 


54+3 


>0.93 


4 13+'-22 


107.6/116 


(26),(27),(28) 


PG1211 + 143 


NLQ 


49.4 


0.171±0.001 


635!«^ 


52+^ 


>0.92 


>4.37 


160.7/139 


(29), (30) 


(*) 6.7 keV K„ 


line is preferred 

















Notes. * Sources for which there is a significant detection of a relativistic ionised 6.7 keV K„ line (see text for details). The first 1 1 sources listed 
in the table belong to the flux-limited sample, d, a, and j8 are the disc inclination angle with respect to the observer, the black hole spin and disc 
emissivity, respectively, of the best-fit parameters of the Kyrline model. Except for the hard X-ray counts where la errors are given, errors in any 
given parameter are given at the 90% confidence level. Also, upper and lower limits in any relevant paremeter are given at the 90% confidence 
level. 

*" SY: Seyfert, Q: quasar, NL: narrow line, BL: broad line, NC: not classified. *2) References where evidence for a relativistic Fe K„, line for a 
particular object is reported (references marked in bold) as well as those where alternative scenarios are provided. The list of references is not 
meant to be complete. 



Refe rences. (1) IDone et al. 20001 (2) [McKeman & Yaqoob 2004] (3) [Steenbrugge et al. 2005[ (4) 'Gondoin et al. 2001 1 (5) IReeves et al.^007 
(6) IBraito et al. 20071 (7) [Wilms et al. (2001)[ (8) [F abian et al. 2002t (9) IVaughan & Fabian 2 004b; ( 10) piller et al. 20Q| 
(11) Ponti et al. 200 61 (12)ptley et al. 2004[ (13)IPounds et al . 2004 , (14) Iwasawa et al. 2004 , (15) Nandra e t al. 19991 (16)|Tumer et al. 20()| 
(17) Tombesi et al. 20071 (18) [Reeves et al. 2004', (19) Markowitz et al. 20 ()9j (20) iGondoin et a l. 20031 (21) iPonti et al. 2009' 



(22)|Page et al. 2003 



(23) IPounds eT al. 2003b; (24) Miller et al. 2007, (25) Vaughan et al. 2004a| (26) Fabian et al. 20091 (27) IFabian et al. 2004 



(28) lGallo et al. 20041 (29) IPounds et al. 2003aL (30) iPounds et al. 2009i 



sistent with this threshold, implying deviations from the standard 
scenario either in terms of very high disc inclination (unlikely 
given the relatively unobscured nature of these AGN), element 
abundances lower than solar or a non-standard accretion disc ge- 
ometry (George & Fabian 1991). 

Last, Table [D.ri in Appendix|D]lists the 20 sources with broad 
line upper limits in the flux-limited sample. In particular, four 
sources (MRK 279, NGC 5548, ESQ 198-G24 and MRK 590) 
exhibit upper limits below 40 eV, arbitrarily chosen to be com- 
parable to the lowest line EW detected. Explaining line EWs be- 
low this value in the framework of standard accretion theories is 
not straightforward, and it requires the use of non-standard sys- 
tem properties. Assuming the above threshold, the non-detection 
fraction of broad lines in the flux-limited sample can then be cal- 
culated as 13% (4/31). Otherwise stated, the presence of a broad 
line cannot be excluded a priori in the remaining 87% of the 
sources in the flux-limited sample (although the line is formally 
detected only in 36% of the sources at a 5cr level). This frac- 
tion is interpreted here as the upper limit to the fraction of AGN 
exhibiting relativistically broadened Fe lines in the parent popu- 
lation. 



4.2.2. Detection fraction and warm absorber 

In some cases, the use of the model component ABSORI could 
be regarded as too simplistic, however necessary for the ap- 
proach considered here. Nevertheless, the influence of the se- 
lected warm absorption component for the baseline model in de- 
termining the detection fraction was tested by running a consis- 



tency check over the 31 sources belonging to the flux-limited 
sample. The test replaces the model component ABSORI by the 
more complex and complete ZXIPCF model and looks for a sig- 
nificant change in the fraction of detected broad lines and the 
average broad line EW. ZXIPCF uses a grid of XSTAR pho- 
toionised absorption models and allows defining a covering frac- 
tion, which has been fixed to I, and column density and ionisa- 
tion parameters of the absorbing material, which have been left 
free during the fitting procedure. The results of this test show 
that the introduction of the ZXIPCF model component does not 
significantly provide a better fit, and yield a detection fraction of 
broad iron K„ lines of 39+15% (12/31). This result is fully con- 
sistent with the 36±14% (11/31) derived when using ABSORI 
for the warm absorption, and concludes that the detection frac- 
tion is not significantly affected by the approach taken in this 
work to model the warm absorber. Out of the sources from the 
flux-limited sample listed in Table |2] only ESO511-G030 does 
not show evidence of a broad line at the 5 sigma level when us- 
ing ZXIPCF. In turn, MRK704 and NGC4593 show a broad line 
at the 5.5cr and 5.3cr levels, respectively, both at the limit of our 
detection threshold of 5 cr. The average line EW derived from 
the detected lines is <EW>=154+112 eV (<EW>,,=73+3 eV), 
again both consistent with the results derived when using the ab- 
sorption model ABSORI (see Section l4.3.1l l. 

4.2.3. Detection fraction and X-ray luminosity 

Eleven out of the 13 lines detected correspond to sources classi- 
fied as Seyfert galaxies, while only 2 come from a source clas- 
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sified as quasar. In a similar manner to the previous section, 
this corresponds to a detection fraction of 13+5% and <6% for 
Seyfert galaxies and quasars respectively. If the exercise is re- 
stricted to sources in the flux-limited sample, the detection frac- 
tions are 36+15% and <92% for Seyfert galaxies and quasars, re- 
spectively. Owing to the low number of lines detected in quasars 
and the low number of quasars present in the flux-limited sam- 
ple, only upper limits to the detection fraction can be derived. 

Figure [3] shows the distribution of the 2-10 keV luminos- 
ity for the FERO sample, highlighting those sources for which a 
significant relativistic broad Fe K„ line detection has been found. 
All broad Fe line detections correspond to objects with luminosi- 
ties below ~ 10'*'* erg s"'. To investigate a possible dependence 
of relativistic broadening on the source luminosity, the follow- 
ing three luminosity bins are defined: LI: Lx < 0.25 xlO'*'* erg 
s-i; L2: 0.25 x lO'*'* < Lx < 1.20 xlO'*'* erg s"*; L3: Lx > 1.20 
X 10"*"* ergs ' . Considering the FERO sample, each bin was cho- 
sen to contain the same number of sources (~50), with no further 
distinctions (see Table[3]). No significant difference is found be- 
tween the number of sources with a relativistic line when com- 
paring the lowest luminosity bin LI to the intermediate lumi- 
nosity bin L2, either using the sources from the FERO (0.7cr) 
or flux-limited (0.5cr) sample. The highest luminosity bin L3 
contains no detections. However, with the uncertainties in the 
detection fraction in the L3 luminosity bin (see Table [3]), it is 
possible to estimate the difference in detection fraction between 
the LI and L3 luminosity bins, at the 2.2cr and 0.9cr levels for 
the FERO and flux-limited samples, respectively. To remove the 
dependency with luminosity and check whether an intrinsic dif- 
ference exist between Seyfert galaxies, and quasars in terms 
of detection fraction, a luminosity range containing the same 
number of Seyfert galaxies, and quasars is selected (L: 0.15 x 
10^*4 < Lx < 2.69x10'*'* erg s"'. Table [3]l. The detection frac- 
tions are 10±7% (4/40) and <10% (2/40) for Seyfert galaxies 
and quasars, respectively, for sources in the FERO sample, and 
18+14% (3/17) and <92% (1/3) for sources in the flux-Hmited 
sample. The low number of detections prevents any statistically 
significant results from being drawn. Table [3] summarizes the 
different detection fraction for the different source and luminos- 
ity class within the whole FERO and flux-limited samples. The 
correlation between the Fe K^ line EW and the hard X-ray lu- 
minosity based on individual source detections and within the 
flux-limited sample is investigated in Section l4~4l This analysis 
is extended to the remainder of the FERO sample in the compan- 
ion paper by Longinotti et al. 



4.3. Average properties of the relativistic broad Fe Ka line 
within the flux-limited sample 

In the next two sections, the average properties of the relativis- 
tic broad Fe Kq, line are discussed, first, in terms of the line 
equivalent width and second, in terms of the disc parameters 
derived from the Kyrline model. To derive meaningful results, 
only the 31 sources in the flux-limited sample are considered; 
therefore, information from the 1 1 sources (the first 1 1 sources 
listed in Table |2]i with a broad Fe K„ line detection and the 20 
sources with an upper limit of the broad Fe K„ line are used (see 
TableED. 



4.3.1. Line equivalent width 

For each source with a significant relativistic broad line detec- 
tion, the measured value of the broad Fe K„ line EW is always 
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Fig. 3. X-ray luminosity distribution in the 2-10 keV energy 
band for the FERO sample. This sample has been split into 
Seyfert galaxies (Sy; 82 sources), quasars (QSO; 67 sources), 
narrow line (NL; 30 sources) and broad line (BL; 60 sources). 
The shaded area corresponds to the luminosity of those sources 
within the flux-limited sample for which a significant relativistic 
broad Fe K„ line detection has been found. The numbers shown 
on the right in each panel are the median values of the corre- 
spondent distribution, also marked by the vertical dashed line. 

Table 3. Fraction of sources in the FERO sample with a sig- 
nificant relativistic broad Fe Kq^ line EW detection for different 
source type and luminosity class. 



Obj. Class/ 
Luminosity 


Sample 
FERO Flux Limited 

(%) (%) 


Sy+QSO 


9+3 (13/149) 


36+14(11/31) 


Sy 
QSO 


13+5(11/82) 
<6 (2/67) 


36+15 (10/28) 
<92 (1/3) 


NL 
BL 


13+9 (4/30) 
12+6 (7/60) 


< 100 (2/4) 
35+18(7/20) 


LI (Sy+QSO) 
L2 (Sy+QSO) 
L3 (Sy+QSO) 


16+7 (8/50) 

10+6 (5/50) 

<2 (0/49) 


44+22 (8/18) 

27+23(3/11) 

<47 (0/2) 


LSy 
LQSO 


10+7 (4/40) 
< 10 (2/40) 


18+14(3/17) 
<92 (1/3) 



Notes. Sy: Seyfert galaxy, QSO: quasar; NL: narrow line, BL: broad 
line; LI: Lx < 0.25 x 10"'* erg s^'; L2: 0.25 x lO""* < Lx < 1.20 x 



10'*" erg s-'; L3: Lx > 1.20 x 10"'* erg s"', where LI, L2 and L3 were 
chosen to contain the same number of sources (~50); L: 0.15 x lO"'* < 
Lx < 2.69 xlO''"* erg s"', L was chosen to contain the same number of 
Seyfert galaxy and quasars (40) (see text for details). 



below 300 eV (see Table |2]l. Considering all the detected lines, 
the mean value of the distribution of broad Fe Kq. line EW is 
143+27 eV, with a \cr standard deviation of 91 eV. It has to be 
kept in mind that the derived value for the mean is likely to over- 
estimate the true mean of the distribution. The explanation can 
be found in the fact that the parent distribution of line EW is 
truncated by the ability to detect weak lines, which can only be 
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Table 4. Weighted mean (<EW>h) and Icr^^ standard deviation, 
and mean (<EW>) and Icr standard deviation, of the Fe K„ line 
EW distribution. 



Source 


<EW>„, 


Cw 


<EW> 


(T 


Class 












(eV) 


(eV) 


(eV) 


(eV) 


SY+QS0(11) 


76+3 


39 


143±27 


91 


SY(IO) 


72±3 


37 


141+30 


95 


QSO(l) 


170±22 





170+22 






Notes. Only sources with an Fe K„ broad line detection within the 
flux-limited sample were used to compute these values. A distinction 
has been made according to source class. The number of sources used 
in each case to derive these values is in parenthesis). 



found in sources with very a high signal-to-noise spectrum (see 
Figure [1]). 

To derive a more representative value of the mean EW, the 
weighted mean of the Fe K^, broad line EW distribution was cal- 
culated using the Icr statistical errors on the line EW for the indi- 
vidual weights. In this way, a weighted mean value of 76±3 eV 
with Icr standard deviation of 39 eV is obtained. This value is 
lower than the mean value derived above, and consistent with it 
to within ~1.3cr. There are two important caveats regarding the 
weighted mean. One, the calculation of the weighted mean as- 
sumes that the individual values of the line EW share the same 
underlying distribution and that such distribution is Gaussian, 
which does not necessarily have to be the case. Second, the use 
of the inverse square of the individual errors in the calculation of 
the weighted mean implies that, for a set of measurements with 
the same relative errors, the weighted mean is biased towards the 
lower values in the set. With a few exceptions, the relative errors 
in the line EW are comparable regardless of the value of the line 
EW, which implies that lower values of the line EW contribute 
more to the weighted mean. As a result, one could consider these 
two values, i.e. the mean and the weighted mean, one biased to- 
wards higher EW values and the other towards lower EW values, 
as the limits where the true mean of the distribution of relativis- 
tic broad Fe K„ line EW in AGN would be expected. Table |4] 
summarizes these results. 

Figure |4] shows the normalised distribution of the line EW 
for sources with a significant broad Fe K^, line detection within 
the flux-limited sample. This distribution is generated by com- 
bining one Gaussian distribution per measurement, where each 
Gaussian takes the mean and standard deviation of the corre- 
sponding best-fit-parameter value and its error at the Icr level 
(upper or lower limits are not present in this plot). Overlaid 
on top of the Gaussian distribution is the corresponding 1- 
dimensional histogram, highlighting the contribution of each in- 
dividual measurement to the overall distribution. Both distribu- 
tions are normalised to their respective peak values. 



4.3.2. Disc parameters 

The Kyrline model parameters 9, /3, and a respectively corre- 
spond to the disc inclination angle with respect to the observer, 
the power-law index of the radial disc emissivity law, and the 
black hole spin, and together they provide important informa- 
tion on the properties and conditions of the accretion disc. These 
values are reported in Table |2] for the sources with a signifi- 
cant broad Fe K^, line detection. For some sources, these pa- 
rameters could not be constrained in the spectral fit and in these 
cases lower or upper limits are reported. Using only the param- 




Fig. 4. Normalised distribution (shaded area) for the line EW for 
sources with a significant broad Fe K^, line detection within the 
flux-limited sample. The 1 -dimensional histograms of this dis- 
tribution is overlaid (white histogram). 



eters from the 1 1 line detections within the flux-limited sample, 
the mean and Icr standard deviation values of 6, /3 and a are 
derived via the censored mean method described in Bianchi et 
al. (2009a). While these authors include measurements and up- 
per limits, lower limits are also considered here. Table |5] lists 
these mean values and the Icr standard deviation obtained using 
this method. (With only two spin measurements, the mean black 
hole spin is omitted.) Figure |5] shows normalised distributions 
for 9 and /?. As for the case of the distribution of the line EW 
in Section 14.3.11 these distributions are generated by combin- 
ing one Gaussian distribution per measurement. The black hole 
spin is constrained well only in two cases, MCG-6-30-15 and 
MRK509, with a best-fit value of 0.86 (!°°^) and 0.78 (!°°^), 
respectively. 

The mean value of /3 is about 2.4, but with a wide spread 
of values. Among the 1 1 sources, 6 have values of /3 higher 
than and inconsistent with 2, which means that emission is 
dominated by the innermost accretion disc regions. In three 
sources /? is definitely lower than 2, and the emission is there- 
fore dominated by the outermost regions. It must be recalled 
that a power law is often a poor approximation of the real emis- 
sivity law, which depends on the (unknown) geometry of the 
emitting corona. (See for instance Martocchia et al. 2000 and 
Martocchia et al. 2002 for the emissivity laws in the lamp- 
post geometry.) It is, however, worth noting that in two sources, 
MCG-6-30-15 and MRK 509, /3 is about 4, or larger, which, un- 
der the lamp-post assumption, would imply an emitting source 
located at only a few gravitational radii from the black hole 
(Martocchia et al. 2002). In addition and for the case of MCG- 
6-30-15, [WihnsetaL(2001) argue that this is not sufficient to 
explain the high yS value and conclude that it is necessary to in- 
voke magnetic extraction of the black hole spin energy to explain 
it piandford & Znajek 1977l l. 

Not much can be said about the spin distribution. In five 
cases this parameter is totally unconstrained, since any permitted 
value of a is acceptable, and in another case it is rather poorly 
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Fig. 5. Normalised distributions (shaded area) for the disc inclination angle (left) and disc emissivity (right) for sources with a 
significant broad Fe Kq, line detection within the flux-limited sample. The 1-dimensional histograms of these distributions are 
overlaid (white histograms). 



Table 5. Mean and \cr standard deviation of the best-fit Kyrline 
model parameters for sources with a relativistic broad Fe Kq, line 
detection. 



Line Parameter 
(Kyrline) 



Mean 



en 

p 



28±1 
2.4±0.1 



5 
0.4 



Notes. With only two spin measurements, the mean black hole spin is 
omitted. 



constrained (a > 0.16). In three cases, lower limits of about 0.5 
can be placed, while in the remaining two sources, the measure- 
ments imply high values of the spin, but inconsistent with maxi- 
mally rotating. It is interesting to note that, whenever a constraint 
can be placed, it always implies the rejection of the static black 
hole solution. However, this is a likely observational bias since 
non rotating black holes produce narrower lines, which are eas- 
ier to detect with respect to the underlying continuum than lines 
produced by maximally spinning black holes, where the red tail 
of these lines can be confused with the underlying continuum. A 
maximally spinning black hole solution is statistically favoured 
than a static black hole solution. Given how few sources there 
are, the approximations made (see for instance the above dis- 
cussions on p) and the subtle interplay between the various disc 
parameters, it is however premature to draw any firm conclusion 
in this respect. 

The relation of the relativistic broad Fe K„ line intensity with 
9, and/? is also plotted in Figure |6] As for black hole spin, con- 
straints are too few (two sources only) to allow for discussion 
of any possible correlation of the line EW with a. These two 
distributions only show two things. First, the distribution of Fe 
K„ line intensity with 9 is consistent with the fact that relativis- 
tic lines arising in accretion discs seen at high inclination angles 
rather than from face-on discs are more difficult to detect for a 
given value of the line EW (as noted by Nandra et al. 2007, see 



also Matt et al. 1992). Second, the apparent lack of detections 
of small EW lines with steep emissivity profiles is most likely 
an observational bias, and any real correlation between the two 
physical quantities can not be claimed. 

4.4. Correlation studies witliin fhe flux-limited sample 

Given all the upper limits derived for the line parameters, the 
censored fit approach was adopted. This procedure for linear 
fitting is applicable to data containing both measurements and 
upper and/or lower limits, and it is based on the method used by 
[Schmitt (1985)| and [Is"obe et al. (1986)| Here it was applied in the 
same way as in Bianchi et al. (2009b). For each correlation anal- 
ysis, the mean Spearman's rank correlation coefficient and null 
hypothesis probability were calculated. Only correlations with a 
null hypothesis probability lower than 10"^ (99.9% c.l.) will be 
considered statistically significant. 

The correlations investigated in this work involve, on one 
hand, the broad Fe K(j, line EW and the disc emissivity (J3), and 
on the other, some source physical properties such as the black 
hole mass (Mbh), the Eddington luminosity (LEdd), and the ac- 
cretion rate (expressed in terms of j^^). Within the censored fit 
method, the errors used for the line EW are Icr errors, while the 
upper limits are at the 90% c.l.. In those correlations involving y6, 
the censored fit approach described in Bianchi et al. (2009b) was 
slightly modified to account for the lower limits to the /3 param- 
eter. The lower limit to /3 can be interpreted simply as consistent 
with 6 at the 90% c.l.. This value of 6 is the maximum value 
allowed for /J in our model (see Section [3.3.1l l. 

The values of the black hole mass were taken from 
the CAIXA catalogue dBianchi et al. 2009al l. and they are 
available for 24/31 of the objects in the flux-limited sam- 
ple. (The reader is deferred to this paper for the references 
from which these values have been extracted.) Using the 
black hole mass (Mbh), the Eddington luminosity was 
derived in this work according to the usual definition: 
Lem = (4nGMmpc) / cr, = 1.26 x 10^** (Mbh/Mq) erg s''. 
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Fig. 6. Dependence of disc inclination angle and emissivity with the relativistic broad Fe K„ line EW. Filled circles indicate line 
detections at >5cr confidence level (where error bars indicate the 90% confidence level intervals), line upper limits are given at the 
90% confidence level. Kyrline parameters upper and lower limits are given at the 90% confidence level. 



The accretion rate is considered as the ratio of the bolometric 
luminosity to the Eddington Luminosity: — —, where the bolo- 

^Edd 

metric luminosity has been derived by applying the following 
bolometric correction: Lboi - khX L^, where k^ is the correction 
factor dependent on the 2-10 keV X-ray luminosity {Lx) follow- 
ing Mar coni et al. (2004)[ The hard X-ray luminosities derived 
in this work for those sources with a broad Fe Kq. line detection 
are reported in Table |6l and they are in good agreement with 
those derived in the CAIXA catalogue. 

Table Q and Figure Q summarize the results of the corre- 
lation studies carried out. Table [T] lists the mean value of the 
Spearman's rank coefficient for the correlation of the broad Fe 
Kq. line EW and (3 with the three physical parameters of inter- 
est, together with the corresponding null hypothesis probability. 
A negative value of Spearman rho indicates anti-correlation be- 
tween the two parameters of interest. On the basis of these re- 
sults, it is concluded that no statistically significant correlation 
is found with any of the considered physical properties of the 
AGNs. 

Last, the relativistic Fe K^, line EW derived in this work and 
all the different parameters included in the CAIXA catalogue 
were correlated (see Table 1 in' Bianchi et al. 2009a| l. No statis- 
tically significant correlation has been found with any of the pa- 
rameters, including the FWHM of the H^ line and the narrow 
Fe 6.4 keV Kq, line EW. However, the low statistics prevents us 
from drawing any robust physical conclusions from the lack of 
such correlations. 



5. Discussion 

Despite the fact that the first relativistically distorted line 
in an AGN was observed about 15 years ago by ASCA 
JTanaka et al. 19951 INandra et al. 199 7). and despite strong ef- 
forts to study them with more recent observatories, like XMM- 
Newton and Suzaku, there are still large uncertainties on the rela- 
tive number of AGNs that possess such lines. The FERO project 



has been designed to address the fundamental question of how 
often these lines appear in AGN by estimating the fraction of 
objects where relativistic effects are detected on a well-defined 
sample. This number is crucial to establish if and how the stan- 
dard AGN paradigm needs to be modified. 

A collection of 149 radio-quiet Type 1 AGN, all targeted 
by XMM-Newton, has been assembled (the largest ever assem- 
bled for this kind of study) in order to tackle the question posed 
above. One fundamental requirement of the FERO project was 
to identify a complete and unbiased subsample of sources with 
high signal-to-noise spectra to be able to derive meaningful con- 
straints on the physical properties of the sources. Only then it 
would be possible to extrapolate the results derived from this 
subsample to the unknown parent population of local radio-quiet 
AGN. A flux-limited sample of 31 objects was defined by se- 
lecting the XSS sources (see Section I2.2l i with count rates in 
the 3-8 keV energy band greater than 1 cts/sec that fulfill the 
FERO selection criteria (i.e. local Type 1 radio-quiet X-ray un- 
obscured sources). In this section the statistical properties of the 
flux-limited sample are discussed. 

The baseline model used in FERO aims at providing the most 
accurate representation of the physical ingredients that are typ- 
ically found in an average Type 1 AGN, so this model includes 
a Compton reflection component and absorption by line-of-sight 
ionised gas for modelling the continuum, a series of narrow Fe 
lines (5), and a relativistic Fe K(j, line component fitted by the 
Kyrline model. The spin of the black hole, the radial emissiv- 
ity profile of the accretion disc, and the disc inclination angle 
to the observer are free parameters in the fit of the line profile. 
To achieve a uniform analysis over the whole FERO sample, no 
attempt has been made to find the best-fit model on a source- 
by-source basis, nor has any attempt been made to find a best-fit 
model by using the FERO sample as a whole. To reduce the pa- 
rameter space during the fitting procedure, several model param- 
eters were investigated to see whether they could be constrained 
due either to physical, geometrical, or observational considera- 
tions. Tests were also carried out on including different model 
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Fig. 7. From top to bottom, relativistic broad Fe K^, line equivalent width andyS vs. black hole mass, accretion rate and hard X-ray 
luminosity. Filled circles indicate broad Fe K„ line detections at >5cr confidence level (where error bars indicate the 90% confidence 
level intervals). The dashed vertical lines in the bottom plots indicate the LI, L2, and L3 luminosity bins as described in the text 
(see Section l4.2.3b . 
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Table 6. Source physical properties used in this work for those sources with a relativistic 6.4 keV Fe Kq, line detection within the 
flux-limited sample. 



Source 




Log,o(BH) 


Lbol 


L'Eddington 


L'bol/^Eddinstoii 


Name 














(10^'' ergs-') 


(Mo) 


(10*^ ergs-') 


(10'"'ergs-') 




IC4329A 


0.56 


8.08 


15.97 


7.42 


0.46 


MCG-5-23-16 


0.14 


7.85 


2.72 


89.20 


0.03 


ESO511-G030 


0.22 





4.87 








MCG-6-30-15 


0.05 


6.19 


0.82 


1.95 


0.42 


NGC4051 


0.003 


6.13 


0.02 


1.70 


0.01 


NGC3516 


0.03 


7.36 


0.49 


28.86 


0.02 


NGC3783 


0.11 


6.94 


2.05 


10.97 


0.19 


NGC3227 


0.01 





0.12 








MRK509 


1.04 


7.86 


35.37 


91.28 


0.39 


MRK766 


0.06 


6.28 


0.93 


2.40 


0.39 


ARK 120 


0.90 


8.27 


29.58 


234.6 


0.13 



Table 7. Results of the correlation studies within the flux-limited sample. 



Broad Fe K„ Line EW (Y) vs. 


a 


b 


Mean Spearman rho 


Null Hypothesis Probability 


Black Hole Mass (X) 
Accretion Rate (X) 
X-ray Luminosity (X) 


1.64±0.12 
1.82±0.15 
7.62±5.16 


-0.20±0.11 
0.08±0.13 
-0.14±0.12 


-0.34 
0.15 
-0.13 


0.14 
0.50 
0.51 


Beta (Y) vs. 


Black Hole Mass (X) 
Accretion Rate (X) 
X-ray Luminosity (X) 


0.36±0.11 
0.38±0.12 
-0.32±4.19 


-0.02±0.09 
0.01±0.10 
0.01±0.10 


0.05 
0.12 
0.17 


0.75 
0.59 
0.39 



Notes. Mean intercept (a) and slope (b) for the linear fits of a log (K) ■ 
are the statistical one standard deviation. 



a + hx log (X) function described in the text. Errors in these two quantities 



components in the baseline model, and the results have served 
either to consider them in the final baseline model or rule them 
out. The most important of these tests was the replacement in 
the baseline model of the neutral reflection component by an 
ionised reflection component, with the consequent change in the 
relativistic line from neutral iron to ionised iron. In fact, the re- 
sults from these two scenarios were merged and presented as if 
coming from a single one. 

Last, a consistency test was carried out by using the hard 
X-ray spectra of all the sources in the flux-limited sample, 
for which a significant detection of the relativistic line can be 
claimed (those listed in Table |2]i, by describing it with the most 
self-consistent reflection model envisaged, however keeping the 
model as simple as possible. The comparison of the results of 
this more complex and self-consistent model for the reflection 
components indicates that it does not significantly affect the re- 
sults derived from the more phenomenological approach chosen 
for the baseline model adopted throughout the paper (refer to 
Appendix|E]for more details). 



5. 1. Fraction of relativistic lines 

The main result of the FERO analysis is the relativistic broad 
line detection fraction, 36 + 14% in the flux-limited FERO sub- 
sample. The observed detection fraction is regarded as a lower 
limit to the intrinsic number of AGN that would show a broad Fe 
K„ line if, for example, all sources were observed with the same, 
and sufficient, signal-to-noise ratio. It is important to stress once 
more that the observed detection fraction comes as a result of 
adopting a rather conservative detection threshold of 5cr con- 
fidence level, which is necessary since FERO does not focus 
on reporting individual broad line detections in particular ob- 



jects. Instead, it focuses on statistical studies of the sample as a 
whole, which includes spectra of very different statistical quality. 
Nevertheless, as a reference, by relaxing this requirement to a 3cr 
confidence level, the above detection fraction rises to 61 ± 20% 
(19/31). 

Upper limits on the line EW can also provide important 
information. For example, the fraction of sources in the flux- 
limited sample in which the broad Fe K„ line EW is limited be- 
low 100 eV is of the order of 50% (where 100% of line EW 
UL are below 360 eV). This result points to a much lower value 
of the average line intensity derived for local Type 1 AGN com- 
pared to the one derived from the mean spectra of high redshift (z 
> 1) Type 1 AGN (<EW> ~ 560 eV; Streblyanska et al. 2005| l, 
although there should be some caution in adopting this result 
since the high mean EW reported could come as a result of 
the method adopted for co-adding the spectra ( Yaqoob2006| l. 
The result derived in this work is also consistent with the up- 
per limit of <400 eV derived by Corral et al. (2008). Using re- 
sults on individual sources, the FERO analysis has provided up- 
per hmits below 40 eV in four sources (MRK 279, NGC 5548, 
ESO 198-G24, and MRK 590), which ai-e constraining enough 
to have to stretch disc-accretion theories by requiring the use of 
non-standard system properties, such as very high disc inclina- 
tion angles or element abundances significantly lower than solar 
( |George& Fabian 199 l| l. 

The non-detection fraction of broad lines in the flux-limited 
sample is 13% (4/31). Inverting this argument, the presence of a 
broad K„ line cannot be excluded a priori in the remainder 87% 
(although the line is formally detected in 36% (11/31) of the 
sources). The value of 87% can be regarded as the upper limit 
to the intrinsic fraction of AGN where a relativistic Fe K„ line 
with EW>40 eV could be present in the parent sample. Nandra 
et al. (2007) have fully explored the range of effects (geometry. 
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ionisation, and element abundances) that would explain why a 
broad line is seen in some objects but not in others, so they are 
not repeated here. 

A direct comparison of the results extracted from FERO 
with previous works reveals that the detection fraction of 36% 
derived for the flux-limited sample is consistent with previ- 
ous results, such as those reported in Nandra et al. (2007) 
where a detection fraction of 35-40% was obtained. When com- 
pared to similar studies on sizable XMM-Newton samples (e.g. 
INandra et al. 20071 Brenneman & Reynolds 2009 1, our sample 
differs in that our sample is very well-defined and homo- 
geneously selected. Moderately X-ray obscured sources and 
radio-loud sources are both excluded from the FERO sample. 
However, the most significant differences are the model used 
to fit the relativistic profile, the fitting procedure, and the treat- 
ment of the sources with multiple observations. For instance, 
Brenneman & Reynolds (2009) fitted the continuum model on 
the entire XMM-Newton energy band, but they selected only one 
particular observation for each source, whereas Nandra et al. 
(2007) include multiple observations of the same target without 
combining them. 



5.2. Disc parameters 

The second goal of the FERO project was to derive average line 
and system properties and look into the question of whether a 
relativistic broad line can be associated to any source physical 
properties and, if so, investigate this dependence in more detail. 
The average broad line EW is of the order of 100 eV, however, 
a wide range of broad line EW (40-260 eV) are observed. As 
a reference, a neutral reflector covering half of the sky, as seen 
from the X-ray source (Q/27r ~ 1) and irradiated by a power- 
law spectrum with F ~ 1.9, should produce an Fe Kq, line with 
equivalent width from around 100 to 150 eV in the incUnation 
range 0°-60° and for solar abundances ( [George & Fabian 1991| l. 
Therefore, line EWs measured outside that range are likely 
to be produced mainly because i) the Fe abundance is differ- 
ent from solar, and/or ii) the reflector (disc) subtends a non- 
standard (smaller) solid angle as seen from the X-ray source. 
One physical scenario where broad line EWs > 150 eV are 
plausible, involves reprocess emission from a photoionised ac- 
cretion disc (' Matt, Fabian & Ross 1993l l or under a geometrical 
scenario which involves anisotropic emission. 

The shape of the broad lines observed in the X-ray spectra of 
AGN are a potential carrier of information concerning the black 
hole spin. The issue of placing constraints to the black hole spin 
has been pursued by several authors (e.g. Nandr a et al. 20071 
Brenneman & Reynolds 2009) l. Using different approaches, both 
studies reached similar conclusions. While some objects show 
evidence of a spinning black hole, with the present data it is very 
difficult to reject the static solution. In this work the approach is 
that of modelling the line with a code capable of including the 
spin as a free parameter The Kyrline model explicitly depends 
on the black hole spin. Kyrline is part of a suit of routines called 
KY (Do vciak et al. 2004] l. a code that properly takes relativistic 
boosting, Ught bending, and time delays into account, assum- 
ing an axysymetric accretion disc geometry surrounding a black 
hole. As a general statement, the black hole spin is poorly con- 
strained by the FERO analysis, except for MCG-6-30-15 (0.86 
!°°i) and MRK 509 (0.78 !°°;^) where the Kerr black hole solu- 
tion is favoured. In five other sources, the static solution cannot 



be rejected, while in a further five, the constraint implies values 
of the black hole spin above ~0.5. 



Another aspect of the present study is to try to establish the 
physical drivers behind the relativistic broad line by looking into 
any correlation between the presence and/or properties of the 
broad line EW and sources physical properties, such as black 
hole mass. X-ray luminosity, or Eddington ratio. However, using 
only sources from the flux-limited sample, no significant cor- 
relation has been observed between the broad Fe K^^ line EW 
and these source properties. As previous works on large samples 
have done (e.g. Iwasawa & Taniguchi 1993 IBianchi et al. 20071 
ICorral et al. 2008l l. the FERO analvsis has looked for inverse 
correlation between Fe Kq. line emission and X-ray luminos- 
ity, the so-called IT effect (previously known a X-ray Baldwin 
effect), where the EW of the narrow Fe Kq, emission line is in- 
versely proportional to the 2-10 keV X-ray luminosity. Earlier 
results on a smaller sample (~100 AGN) using the spectra 
stacked residuals technique suggested that relativistic broad pro- 
files follow a similar trend, where broad profiles are more com- 
mon in low-luminosity AGN JGuainazzi et al. 20061) . With infor- 
mation from individual source line profiles, the FERO analysis 
is unable to confirm such results, finding no statistically signif- 
icant indication for the IT effect on broad Fe lines. However, it 
is worth pointing out that the lack of sources in the flux-limited 
sample with luminosities above lO**"* erg s"' prevents any firm 
conclusions. The IT effect will be investigated using the spec- 
tra stacked residuals technique using those sources in the FERO 
sample with an upper limit to the relativistic Fe Kq. line EW 
(Longinotti et al. in preparation). 



5.3. FERO in a wider context 

An important remark should be made about the relativistic line 
in sources that are characterized by other phenomena. A recent 
paper by |De Marco et al. (2009)| carried out a systematic search 
for the presence of red- or blue-shifted components in the time- 
averaged spectra of the sources in the flux-limited sample. The 
study also looks into the spectral variability of these spectral 
components. A direct comparison with their results shows that 
10 out of the 1 1 sources (~90%) within the flux-limited sample 
(all except ESO 511-G030) with a relativistic broad line indi- 
cate either red- or blue-shifted components, or both. In 6 out of 
these 10 sources (IC4329A, NGC 3783, NGC 3516, MRK 509, 
MCG-6-30-15, and MRK 766) there is evidence of significant 
variability (>90%) in the excess map (5.4-7.2 keV). 

From the literature on the warm absorber phenomenon, 
we note that 9 out of 11 sources with line detections have 
warm absorbers: NGC3783, N GC3516, MRK 50 9, MCG-6- 
30-15, MRK 766, NGC 4051 ([Blustin et al. 20051 McKernan 
et al. 20 07), IC4329A (jSteenbrugge et al. 2005|l, NGC 3227 
(Markow itz et al. 20091) . In the case of ESO 5 1 1-G030 the pres- 
ence of line-of-sight ionised absorption has not been ascertained 
yet (Longinotti private communication). The only 2 sources in 
the flux-limited sample with no warm absorber and no transient 
features are Ark 120 and MCG-5-23-16, so it is tempting to 
deem these two sources as the most robust detections, due to 
the lack of an alternative scenario to the relativistic line. 

That sources with significant relativistic lines simultaneously 
show Fe K band variability and a warm absorber is probably 
twofold. There might be a physical connection between, e.g., 
the presence of the relativistic line and of the transient features, 
with some variability associated to it. But, it could well be the 
result of an observational bias, because these sources were ex- 
tensively (and repeatedly) observed because they were known to 
be interesting a priori, so they ended up with the best signal-to- 
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noise data and this explains the coincidence of the broad Fe line, 
warm absorber, and transient features in the same sources. 

To conclude, it goes without saying that the FERO project 
would benefit from an expansion in several directions. The most 
obvious would be to complete the observations of the flux- 
limited sources such that the number of X-ray counts in the 2- 
10 keV fall above 1.5x10^ cts. The expansion of the number of 
sources with 2-10 keV X-ray luminosities above lO'*'* erg s"' 
would also allow studies like the IT effect to be carried out at a 
significant level and, of course, to investigate other correlations 
where the X-ray luminosity plays an important role. 

Another crucial aspect to consider in future studies regard- 
ing relativistic disc lines is broadband coverage. This would al- 
low exploration of the spectral regime above 10 keV and the op- 
portunity to discern whether the spectral contribution at these 
energies can be attributed to a reflection component off dis- 
tant matter and/or the accretion disc, analogous to the case of 
MCG-5-23-16 JReevesetal. 20071 ). The importance of broad- 
band coverage is also highUghted by, for example, the recent 
detection of a broad Fe K^, line in Fairall 9 in Suzaku observa- 
tions, which have allowed a tight constraint on the black hole 
spin dSchmoU et al. 20091 1. 

The FERO project has highlighted the importance of the 
availability of a complete and well-exposed AGN sample to shed 
some light on the origin and properties of the relativistic broad 
Fe lines in AGN, the environment, and conditions in which they 
are produced, and to provide better understanding of the scien- 
tific drivers needed for future high-throughput X-ray missions, 
such as IXO. 



6. Conclusions 

The main conclusions that are drawn from the results of the sys- 
tematic and homogeneous spectral analysis of the large number 
of sources in FERO can be summarised as follows. 

• The observed fraction of sources in the FERO sample that 
present strong evidence (>5cr significance) of a relativistic 
broad Fe K„ line is of the order of 9% (13/149). Considering 
only sources from the flux-limited sample, the detection frac- 
tion rises to around 36% (11/31). This number can be in- 
terpreted as a lower limit to the intrinsic fraction of AGN 
with relativistic Fe line broadening. The sources with a sig- 
nificant Fe K„ line detection mostly belong to sources with 
good statistical quality spectra (> 1.5x10^ cts), so it is not 
surprising that the majority belong to the flux-limited sam- 
ple. Considering only sources from the flux-limited sample, 
a broad Fe Kq, line at the level of 40 e V can be rejected in 4 
objects, 87%, which can be regarded as an upper limit to the 
intrinsic fraction of AGN with relativistic Fe Kq, line broad- 
ening. 

• There is no significant difference between Seyferts and 
quasars or narrow line and broad line objects in terms of de- 
tection fraction. 

• All sources with a broad Fe K^, line detection have a hard X- 
ray luminosity in the 2-10 ke V energy band below ~ 1 x 1 0"*"* 
erg s"^ The difference in the detection fraction between the 
highest and lowest luminosity bins as defined in this work is 
only 2cr. 

• From those sources with a significant Fe K„ line detection 
within the flux-limited sample, it is found that: 

o The average relativistic Fe Kq, line EW is of the order of 
100 eV (never higher than 300 eV for any given source). 



o The average system properties and Icr standard deviation 
inferred from the Kyrline model can be summarised as: 

< 6* > = 28 ± 5°, consistent with an intrinsic random 
distribution of inclination angles. 

< /3 > - 2.4 ± 0.4, with a wide spread of values. 

The spin value a is poorly constrained, except for MCG- 
6-30-15 (0.86 (+;;;;;')) and MRK509 (0.78 (!°;°^)) where 
the Kerr black hole solution is favoured. 

o No significant trend is found between the Fe K„ line EW 
and any of the line parameters investigated (9, /3, and a). 

o No significant correlation has been found between the Fe 
Kq line EW or disc emissivity (j3) and the source physical 
properties investigated, such as, black hole mass, accre- 
tion rate, and hard X-ray luminosity. 
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Table B.l. List of sources within the FERO sample where mul- 
tiple observations are available but only one has been used (the 
one with longest exposure time). 



Appendix C: Best-fit-model parameters for sources 
belonging to the flux-limited sample. 



Source 


List of 


Date 


Exposure 




XMM-Newton 




(ksec) 




Observation 








IDs combined 






1H0707-495 


0148010301 


2002-10-13 


70.0 




0110890201 


2000-10-21 


36.4 


ESO113-G010 


0301890101 


2005-11-10 


90.6 




0103861601 


2001-05-03 


4.0 


ES0198-G24 


0305370101 


2006-02-04 


85.3 




0067190101 


2001-01-24 


26.4 




0112910101 


2000-12-01 


6.8 


IZWl 


0300470101 


2005-07-18 


58.0 




0110890301 


2002-06-22 


18.4 


MR225 1-178 


0012940101 


2002-05-18 


44.5 




0112910301 


2000-11-29 


3.5 


MRK841 


0205340201 


2005-01-16 


30.1 




0112910201 


2001-01-13 


5.9 




0070740101 


2001-01-13 


7.6 




0070740301 


2001-01-14 


9.3 




0205340401 


2005-07-17 


17.3 


NGC2622 


0302260201 


2005-04-09 


6.6 




0302260401 


2005-10-08 


4.5 


NGC4051 


0109141401 


2001-05-16 


73.3 




0157560101 


2002-11-22 


46.2 


NGC4593 


0059830101 


2002-06-23 


53.2 




0109970101 


2000-07-02 


8.8 


PG1211 + 143 


0112610101 


2001-06-15 


49.5 




0208020101 


2004-06-21 


43.6 


UGC3973 


0400070201 


2006-09-30 


14.5 




0103860801 


2000-10-09 


1.7 




0103862101 


2001-04-26 


3.6 


PG1116+215 


0201940101 


2004-12-17 


69.3 




0201940201 


2004-12-19 


5.0 




0111290401 


2001-12-02 


5.5 


IRASF12397+3333 


0202180201 


2005-06-20 


68.7 




0202180301 


2005-06-23 


9.1 



Notes. Observation IDs in bold correspond to those observations that 
have been considered, while observation IDs not in bold have been dis- 
carded. 



Appendix A: Sources from the flux-limited sample. 

List of 33 sources in the flux-limited sample. These sources 
have been selected from the RXTE all-sky Slew Survey (XSS, 
jRevnivtsev et al. 2004) with a count rate in the 3-8 keV energy 
band greater than 1 cts/sec and fulfilling the FERO source selec- 
tion criteria. The XSS is nearly 80% complete at the selected flux 
level for sources with Galactic latitude greater than 10°. For two 
of them, UGC 10683 and ESQ 0141-G055, no XMM-Newton 
data were available as of April 2008. This leaves the number of 
XSS-selected bright sources in the FERO sample at 3 1 . 



Appendix B: List of sources with multiple 
observations that have not been summed. 

List of sources within the FERO sample where multiple obser- 
vations are available but only one has been used (the one with 
longest exposure time). 



Table IC.ll gives some relevant best-fit-model parameters (see 
Section [33T l for the 3 1 sources belonging to the flux-limited sam- 
ple. 



Appendix D: Relativistic Fe K,, line EW upper limits 
for the sources in the flux-limited sample. 

List of 20 sources within the flux-limited sample with an upper 
limit to the relativistic Fe K^ line EW. 



Appendix E: A final self-consistent test on the 
detections within the flux-limited sample 

As a final test, the hard X-ray spectra of all the sources in 
the flux-limited sample for which a significant detection of 
the relativistic line can be claimed has been described with 
the most self-consistent reflection model envisaged, but keep- 
ing the model as simple as possible. The baseline model com- 
prises one layer of ionised absorption (the ZXIPCF model), re- 
flection off cold distant matter including the most important 
emission lines and the associated self-consistent reflection con- 
tinuum (the PEXMON model), and reflection from the accre- 
tion disc (the IRoss & FabianI REFLION model, convolved with 
the KYRLINE kernel). Two ionised emission lines with energies 
fixed at 6.7 keV and 6.96 keV are also included as in the previous 
phenomenological models used throughout the paper, although 
they are not statistically required in all cases. 

The PEXMON model is described in detail in Nandra et al. 
(2007) and describes the reflection spectrum from a cold slab 
of gas including both the reflection continuum and the most rel- 
evant emission lines, which are computed as self-consistently 
as possible according to the work by George & Fabian (1991). 
The metal abundances are fixed to the solar value (except for 
the case of MCG-6-30-15, see below) and the reflector incli- 
nation with respect to the line of sight to 60 degrees, which is 
appropriate for torus-like reflection in Seyfert 1 galaxies. The 
illuminating continuum is a power law with the same photon 
index as the power-law component of our spectral model. The 
REFLION model describes reflection off a ionised slab and is 
used instead to describe the disc reflection component. The illu- 
minating photon index is the same as the primary continuum and 
the Fe abundance is fixed to the solar value (except for MCG-6- 
30-15, see below). The model is convolved with the kernel of the 
KYRLINE model (Dovciak et al. 2004), which allows including 
all relativistic effects and measuring the relevant parameters. 

It is important to stress once again that the goal here is not 
to provide the best possible fitting statistics but rather to com- 
pare the results of the phenomenological model used in the paper 
(which does not account for emission lines and associated reflec- 
tion continua in a self-consistent manner) with a more physical 
spectral model. In some cases, more spectral components are in- 
cluded as explained in the subsequent section. 

The results are reported in Table IE. ll where the most impor- 
tant spectral parameters associated with the reflection compo- 
nents are considered. In Table IE. 2l the best-fitting relativistic pa- 
rameters for both the phenomenological and the more physically 
motivated models are reported. Such a comparison implies that 
the more complex and self-consistent model for the reflection 
components does not significantly affect the results. Since the 
cases considered here correspond to the highest signal-to-noise 
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Table A.l. Flux-limited sample: 33 sources from the RXTE Slew Survey with a count rate in the 3-8 keV energy range greater than 
1 cts/sec and fulfilling the FERO source selection criteria. 



Source 


Type 


RXTE Slew Survey 


XMM-Newton 


XMM-Newton Exposure Time 






(cts/sec)3-8 keV 


(cts/sec)2-io keV 


(ksec) 


NGC4593 


BLSY 


1.05 ±0.16 


4.613 ± 0.009 


53.2 


MRK704 


BLSY 


1.06 ±0.22 


0.984 ± 0.009 


14.9 


ESO511-G030 


NCSY 


1.10 ±0.09 


2.284 ± 0.005 


76.2 


NGC7213 


BLSY 


1.12 ±0.16 


2.448 ± 0.009 


29.6 


AKN564 


NLSY 


1.13 ±0.02 


2.246 ± 0.005 


84.1 


HI 846-786 


NCSY 


1.15 ±0.09 


0.700 ±0.011 


6.3 


MRKUO 


NLSY 


1.17 ±0.10 


3.362 ± 0.010 


32.8 


ESO198-G024 


BLSY 


1.21 ±0.19 


1.150 ±0.003 


85.3 


FAIRALL9 


BLQ 


1.25 ± 0.02 


1.313 ±0.007 


25.9 


UGC3973 


BLSY 


1.32 ±0.12 


1.304 ±0.019 


3.6 


NGC4051 


NLSY 


1.49 ±0.02 


2.783 ± 0.007 


73.3 


NGC526A 


NCSY 


1.61 ±0.14 


2.177 ±0.007 


41.5 


MCG-2-58-22 


BLSY 


1.61 ± 0.08 


3.459 ± 0.021 


7.2 


NGC7469 


BLSY 


1.69 ±0.01 


3.343 ± 0.005 


142.6 


MRK766 


NLSY 


1.77 ±0.12 


1.962 ±0.002 


398.3 


MRK590 


BLSY 


1.95 ±0.10 


0.766 ± 0.003 


72.0 


IRAS05078+1626 


NCSY 


2.08 ±0.19 


2.691 ± 0.008 


40.1 


NGC3227 


BLSY 


2.10 ±0.02 


3.635 ± 0.006 


97.4 


MR225 1-178 


BLQ 


2.10 ±0.11 


2.387 ± 0.007 


44.5 


MRK279 


BLSY 


2.13 ±0.06 


3.110 ±0.006 


104.8 


ARK120 


BLSY 


2.14 ±0.03 


4.553 ± 0.008 


76.3 


NGC7314 


NCSY 


2.16 ±0.04 


4.509 ± 0.012 


30.1 


H0557-385 


NCSY 


2.33 ± 0.34 


4.254 ± 0.021 


10.5 


MCG+8-11-11 


BLSY 


2.36 ± 0.27 


4.838 ±0.014 


26.5 


MCG-6-30-15 


BLSY 


3.08 ± 0.02 


4.719 ±0.004 


295.1 


MRK509 


BLQ 


3.12 ±0.10 


4.432 ± 0.005 


27.8 


NGC3516 


BLSY 


3.21 ±0.01 


1.841 ±0.004 


153.0 


NGC5548 


BLSY 


3.58 ± 0.02 


4.191 ±0.007 


74.8 


NGC3783 


BLSY 


4.90 ± 0.04 


5.539 ± 0.006 


170.4 


MCG-5-23-16 


NCSY 


6.46 ±0.12 


8.025 ± 0.008 


110.3 


IC4329A 


BLSY 


7.29 ± 0.07 


10.713 ±0.011 


85.1 


UGC10683 


SY 


1.27 ±0.18 








ESO141-G055 


SY 


1.37 ±0.12 









Notes. These sources, with the exception of the last two sources, have all been the target of XMM-Newton observations. The XMM-Newton EPIC- 
pn count rate in the 2-10 keV energy range and exposure time is also given. For two of them, UGC 10683 and ESO 0141-G055, no XMM-Newton 
data were available as of April 2008. This reduces the number of XSS-selected bright sources in the FERO sample to 31. 



data within the whole sample, the test supports the analysis car- 
ried out on the whole available sample within the context of the 
less sophisticated and more phenomenological model. 



E.1. Notes on individual sources 

Slight modifications to the baseline model described in 
Appendix |E] are given here for some individual sources. In the 
remaining cases, the baseline model was applied without being 
modified. 

IC4329A: Absorption is best modelled with a neutral layer 
with moderate column density (^ 4 x 10^' cm"^). A Gaussian 
absorption line with EW ^ -15 eV was also included to model 
an absorption feature at ^ 7.65 keV phenomenologically. 

MCG-5-23-16: Absorption is best modelled with a neutral 



layer with column density Nh ^ 1.5 x 10 cm . The presence 
of a further ionised layer is possible but not statistically required. 

MCG-6-30-15: The Fe abundance of the two reflectors 
was left free to vary in this case resulting in an overabundance 
of ~ X 3 with respect to the solar value. The model also com- 



prises a Gaussian absorption line at ^ 6.7 keV with EW ^ -15 eV. 

NGC4051: A Gaussian absorption line with EW^ -30 eV 
was also included to model an absorption feature at ^ 7.04 keV 
phenomenologically. 

NGC3516: A Gaussian absorption line with EW^ -28 eV 
was also included to model an absorption feature at ^ 6.7 keV 
phenomenologically. 

NGC3783: A Gaussian absorption Une with EW^^ -20 eV 
was also included to model an absorption feature at - 6.7 keV 
phenomenologically. 

Mrk509: A Gaussian absorption line with EW- -10 eV 
was also included to model an absorption feature at ^ 7.3 keV 
phenomenologically. 
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Table C.l. Flux-limited sample: relevant best-fit-model parameters corresponding to the neutral reflection run and 6.4 keV relativis- 
tic Fe Kq, line. 



Source 


Power Law 


Warm Absorber 


Compton 


Reflection 


EW Narrow Fe K^ Emission Lines 




r 


Nh 


f 


R 


Norm 


Fe, 6.4 keV 


Fexxv 6.7 keV 


Fexxv/ 6.96 keV 






(10^2 cm~2) 


(erg cm ' s ' ) 


(10- 


3 keV-' cm-2 s"') 


(eV) 


(eV) 


(eV) 


AKN564 
NGC526A 


2 82+"-"^ 
^■■'^-0.02 


76+"''' 

124+8:^^ 

0.80!|?] 


<0.47 


9 99+0.62 

''■^ -1.25 

<0.77 


18.38!^JS 
4.46!iJ;^? 


< 18 

49+13 
4Z_,7 


13+10 
"-11 
<26 


<4 

<22 


MCG-6-30-15 


2 n'i+"-"2 

^■"^-0.01 


6 30+1 -'7 


n<:,+0.29 
U.31 17 


17.38!j>:g 


44+^ 


<5 


<7 


MRK766 


2 iQ+0.06 

1.64+llf 


o.40!j;-|« 


10.37^8-Vi 


o.78!!;-]^ 


>6.87 


32- 


<21 


20!?, 


NGC5548 


<0.42 





0.20+2-^2 


8.57+j;-« 


67+], 


<4 


<7 


MCG+8-11-11 


1 70+0:05 

1-^^-0.04 


<6.24 





, o7+o:M 

' -''-0.48 


11.78!^ 
3.37!ijf3 


97+10 
^'-11 


< 10 


< 17 


MRK704 


1 72+"-"'* 


9.67+^1 
4.48!: 

'■20 002 


ID'S 49+74.83 
1U3.4Z 


<3.33 


<83 


<42 


<45 


NGC3516 
MCG-5-23-16 


1.61+""* 
1.65+llf 


77 03+li-«8 
<0.02 


2 99+0.22 
1.15+IJ^ 


4.02+"-!^ 
19.37"+of« 


140!' 

32+ 


<7 
<5 


<5 
<5 


NGC3227 
H0557-385 


1 1:9 +0:01 
1.86 Qgj 


0.09+"-"' 
1 99+0:57 

^■^ -0.32 


<25.38 
5.67!-- 


0.15+i 
<3.28 


6 70+lfof 

-0,19, 

13.63!*:[2 


59+"' 
<31 


<5 

< 14 


8+1 
<26 


NGC4593 


1 84+""^ 
'■°^-0.04 


<4.07 




1 24+"-32 


ii.88!S:« 


39+26 


< 10 


14!^ 


NGC7469 


>1.97 


<0.53 





1 79+0.25 


10 07+0.20 


73+5 


<3 


<9 


NGC3783 


1 70+"" 1 


1 64+""^ 
'■°^-0.06 


54.21!^";; 


51+""" 
U.31 Q^3 


1 44+0.08 
'^■^^-0.17 


90!^ 


< 1 


23!^ 


ESO511-G030 


1 Qq+0.02 
i.yj 03 


<1.38 





2 in+"-''2 


6.31+0.26 

5 96+8:?^ 


11!!; 


<14 


< 12 


IRAS05078+1626 


1 79+0.03 
'■'^-0.06 


4 94+3.00 
^■^^-4.05 





1 7n+"-28 


38+21 


11+'" 


< 11 


MRK279 
MRK509 


1 Qq+0.02 
1.83_oo3 

>1.96 


<0.51 
<0.2 





1.16+"'" 

2 49+8:35 
■^■^^-0-36 


7 79+0.07 

12.88+j;-|^ 


,^+16 
^^-14 

34!^ 


<3 
<4 


<3 


NGC7314 


2 09+0-11 


1 3S+"-32 


4 86+5-'*8 
^■™-3.8l 


1 16+'"" 
'■'"-0-64 


18.22!J| 


<56 


14+'" 
-!9 


27+11 


NGC7213 


<7.67 





<0.77 


5.80!iJJ« 


82!]^ 


92+13 
zz_,3 


17!|3 


IC4329A 

MR2251-178 
H1846-786 


J -79+0.02 
1 09+0:03 


f) 3^+0.04 
"■^-'-0.02 

<0.61 
<8.95 


<0.47 


0.93+"-°' 

54+":2* 
"■-'^-0.48 

>0 


26 85+"-'* 
306+"""''"' 

: 11 


34+^ 
23+11 

78+6" 


<4 
<5 

<35 


8!^ 

<2 

104+]' 


FAIRALL9 


1.90+™^ 


<2.67 




1.70+?-« 


3.56+f 


107+1 


<21 


-72 

<29 


UGC3973 


1 78+8: 


<37.79 





-1.05 
>0 


3 06+'-^* 
7:65+8:5^ 
9.48+^f 

'■™-0.04 


156!^ 


<87 


<38 


MCG-2-58-22 


1.69+"'" 
202+8:8^ 

1 61+8:8i 
'■"^-0.02 


<2.93 





<3.32 


<51 


<20 


21+27 
< 17 
41!- 


NGC4051 

MRK590 


2 21+3-39 
^■^ -0.64 

<3.39 





<1.49 

59+"-''2 
"■-'^-0.39 


63!|2 
113!!^ 


17+'" 
"-11 

<31 


ARK120 


9 19+0.02 


<0.52 





2 14+0.35 

^■^^-0.34 


15.95!«ii 


65!] 


15!^ 


9!] 


MRKllO 


1 Qs+0.02 


<1.21 





1 47+0.25 
^■^'-0.61 


8.58!iJ;^2 


30+10 


<6 


<3 


ES0198-G24 


1 7Q+0.07 
'■'^-0.04 


<2.57 





, 70+0.83 


9 7+0.16 
^■'-0.09 


60+ '2 


<8 


<8 



Notes. Errors are given at the 90% confidence level. Upper limits to the relevant parameters are also given at the 90% confidence level. Sources 
with a significant broad Fe Kq. line detection are marked in bold (as listed in Table[2ll. 

Table D.l. Flux-limited sample: list of 20 sources with an upper limit to the relativistic Fe Kq, line EW. 



Source 


Cts^-i"^^^ 


EW Upper Limit 




(10^ cts) 


(eV) 


MRK704 


0.146 ±0.012 


362.1 


FAIRALL9 


0.341 ± 0.018 


320.4 


AKN564 


1.889 ±0.043 


292.3 


MCG-2-58-22 


0.250 ±0.016 


208.5 


HI 846-786 


0.044 ± 0.007 


173.7 


MCG-f8-ll-ll 


1.284 ±0.036 


158.4 


NGC526A 


0.903 ± 0.030 


122.1 


MRKllO 


1.104 ±0.033 


115.6 


NGC7213 


0.724 ± 0.027 


107.4 


UGC3973 


0.389 ± 0.020 


100.2 


H0557-385 


0.446 ± 0.021 


99.3 


MR2251-178 


1.062 ± 0.032 


99.2 


NGC7314 


1.358 ±0.037 


87.8 


NGC4593 


2.455 ± 0.049 


84.2 


NGC7469 


4.769 ± 0.069 


84.0 


IRAS05078-t-1626 


1.079 ± 0.033 


61.4 


MRK279 


3.259 ± 0.057 


38.2 


NGC5548 


3.133 ±0.056 


37.6 


ES0198-G24 


0.981 ±0.031 


35.8 


MRK590 


0.552 ± 0.023 


33.4 



Notes. The EW upper limit is provided at the 90% confidence level. 



21 



I. de la Calle Perez et al.: Statistics of Relativistic Fe K„ lines in Type 1 AGN 



Table E.l. Summary of the results of the self-consistent reflection model. 



Source 



Distant reflection Disc reflection 



Disc ionisation 
(erg cm s"') 



IC4329A 




0.5 ± 0.2 


MCG-5-23- 


16 


0.3 ±0.1 


ESO511-G030 


0.7 ± 0.2 


MCG-6-30-15 


0.2 ±0.1 


NGC4051 




0.9 ± 0.3 


NGC3516 

NGC3783 

NGC3227 

MRK509 

MRK766 

ARK120 




1.5 ±0.3 
0.8 ± 0.2 
0.6 ± 0.2 
0.3 ±0.1 
0.3 ±0.1 
1.0 ±0.3 



0.3: 
0.5: 
0.8: 
2.4: 
1.0: 
1.8: 
0.8: 
0.5: 
1.0: 
0.8: 
0.7: 



0.2 
0.3 
0.4 
0.8 
0.5 
0.5 
0.4 
0.3 
0.6 
0.4 
0.4 



<70 
70 ±50 

62^'*'' 

<70 
<90 

-7Q+60 

80!| 
700!^^° 



6 

n 



xVdof 



34+11 

31+^ 
37+21 

39 ±5 

18!' 
34!' 

97+19 
^'-20 

26 ±9 

>45 

28 ±8 

45 ± 10 



>0 

>0 

>0 

>0.8 

>0.3 

>0.3 

>0.2 

>0 

>0.4 

>0.3 

>0.1 



T 4+0.8 
-'•^-0.5 
9 c+0.6 
^•-^-0.5 

2.7±0.4 
3.1!"^ 
2.3+!-! 
2.8 
2.6; 
2 7 



197.2/163 
247.3/164 
185.0/162 
192.1/164 
172.0/162 
212.7/161 
221.4/163 
226.5/162 
244.2/160 
262.5/162 
188.7/160 



Notes. All sources for which a detection of a relativistic line component is claimed in the paper have been described with two self-consistent 
reflection models, one for distant reflection (the PEXMON model), the other for disc reflection (the REFLION model). The relevant parameters 
and fitting statistics are reported in this table. 

Table E.2. Comparison of the phenomenological and self-consistent model. 



Source 



IC4329A 

MCG-5-23-16 

ESO511-G030 

MCG-6-30-15 

NGC4051 

NGC3516 

NGC3783 

NGC3227 

MRK509 

MRK766 

ARK 120 



Phenomenological Model 



9 

n 



/3 



28!T" 

21+s 

18 ±7 
40!J 

22 ±6 
27+2 

<8 

23 ±4 
53 ±1 

20!' 
>59 



>0 


< 1.3 


>0 


< 1.6 


>0 


< 1.1 


OMtZ 


4.1±0.2 


>0.46 


9 Q+0.3 
^•='-0.4 


> 0.48 


2.8!- 


>0.16 


2.7!- 


>0 


1 Q+0.6 
^•='-0.5 


0.78!»- 


> 3.8 



>0.47 
>0 



9 7+0.2 
^•'-0.1 
9 9+0.6 
■^- -0.3 



Self-consistent Model 



6 

n 



p 



31!'' 
37+i? 

J/_20 

39 ±5 
18!?3 
34!' 

97+15 
^'-20 

26 ±9 

>45 

28 ±8 

48 ±12 



>0 

>0 

>0 

>0.8 

>0.3 

>0.3 

>0.2 

>0 

>0.4 

>0.3 

>0.1 



TDT^ 



1.9 



2.6!°-^ 
3 4+Ai 

2.7±0.4 

-'•^-0.4 

2 3+'-' 
r,-l 

2 g+8:I 

-Q.5 
2 7+0.5 



Notes. This table shows a comparison of the best-fitting relativistic parameters (only disc parameters are shown) obtained by applying the 
phenomenological model discussed throughout the paper (whose full results are given in Table |2ll and the more self-consistent one discussed in 
Appendix|E] Here, only sources for which a detection of a relativistic line component is claimed in the paper are shown. 
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